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Abstract 

We study the critical crossover between the Gaussian and the Wilson-Fisher fixed point 
for general 0(AQ-invariant spin models with medium-range interactions. We perform a 
systematic expansion around the mean-field solution, obtaining the universal crossover 
curves and their leading corrections. In particular we show that, in three dimensions, the 
leading correction scales as R~ 3 , R being the range of the interactions. We compare our 
results with the existing numerical ones obtained by Monte Carlo simulations and present 
a critical discussion of other approaches. 



1 Introduction 



Every physical situation of experimental relevance has at least two scales: one scale is intrinsic to 
the system, while the second one is related to experimental conditions. In statistical mechanics 
the correlation length £ is related to experimental conditions (it depends on the temperature), 
while the interaction length (Ginzburg parameter) is intrinsic. The opposite is true in quantum 
field theory: here the correlation length (inverse mass gap) is intrinsic, while the interaction 
scale (inverse momentum) depends on the experiment. Physical predictions are functions of 
ratios of these two scales and describe the crossover from the correlation-dominated (£/G or 
p/m large) to the interaction-dominated or p/m small) regime. In a properly defined limit 
they are universal and define the unique flow between two different fixed points. This universal 
limit is obtained when two scales become very large with respect to any other (microscopic) 
scale. Their ratio becomes the (universal) control parameter of the system, whose transition 
from to oo describes the critical crossover. 

In this paper we will consider the crossover between the Gaussian fixed point where mean- 
field predictions hold (interaction-dominated regime) to the standard Wilson-Fisher fixed point 
(correlation-dominated regime). In recent years a lot of work has been devoted to understanding 
this crossover, either experimentally [|l|-|5j or theoretically |6|-f2"0|j. The traditional approach to 
the crossover between the Gaussian and the Wilson- Fisher fixed point starts from the standard 
Landau-Ginzburg Hamiltonian. On a d- dimensional lattice, it can be written as 



+ E 



U 

\ r 4>l + ^4>t ~h x -(j)a 



'l.F 



where <fi x are iV-dimensional vectors, and J{x) is the standard nearest-neighbour coupling. For 
this model the interaction scale is controlled by the coupling u and the relevant parameters are 



the (thermal) Ginzburg number G |2T| and its magnetic counterpart |T3,^3 defined by 



G 



u 



2/(4-d) 
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(d+2)/[2(4-d)] 



1.2) 



Under a renormalization-group (RG) transformation G scales like the (reduced) temperature, 
while Gh scales as the magnetic field. For t = r — r c ^ G and h ^ Gh one observes the standard 
critical behaviour, while in the opposite case the behaviour is classical. The critical crossover 
limit corresponds to considering t, h, u — > keeping t = t/G and h = h/Gh fixed. This limit is 
universal, i.e. independent of the detailed structure of the model: any Hamiltonian of the form 
Ql.lj) shows the same universal behaviour as long as the interaction is short-ranged, i.e. for any 
J(x) such that X^x 2 J(x) < +oo. The crossover functions can be related to the RG functions 
of the standard continuum 4 theory if one expresses them in terms of the zero-momentum 
four-point renormalized coupling g ]^-^,^T|. For the observables that are traditionally studied 
in statistical mechanics, for instance the susceptibility or the correlation length, the crossover 
functions can be computed to high precision in the fixed- dimension expansion in d = 3 



Let us now consider the medium-range case. Following Refs. [|T^,[T^] we assume that J(x) 
has the following form 
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where D is a lattice domain characterized by some scale R. Explicitly we define R and the 
corresponding domain volume Vr by 



VJ^£i, ^ 2 -^rE- 2 - (1-4) 

The shape of D is irrelevant for our purposes as long as Vr ~ R d for R — > oo. The constant J 
defines the normalization of the fields. Here we assume J = 1/Vr, since this choice simplifies the 
discussion of the limit R —>■ oo. To understand the connection between the theory with medium- 
range interactions and the short-range model let us consider the continuum hamiltonian that 
is obtained replacing in Eq. ( |1 . 1|) the lattice sums with the corresponding integrals. Then let 



us perform a scale transformation [p^| . We define new ("blocked") coordinates y = x/R and 
rescale the fields according to 

= R d/2 Ct> Ryi hy = R^hRy. (1.5) 

The rescaled Hamiltonian becomes 



d d yi d d y 2 \J{ Vl - y 2 ) ($ yi - $ y2 ) + J d d y ^r$ 2 y + -—$ y -h y 
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where now the coupling J(x) is of short-range type in the limit R —>■ oo. Being short-ranged, 
we can apply the previous arguments and define Ginzburg parameters: 

G = {uR- d f {d ~' ] = ^-4)^/(4-^ (L7) 

G h = R- d ' 2 [uR- d ) [d+2)md ~ i)] = u^ d+2 ^ d -^ R-^- d K (1.8) 

Therefore, in the medium-range model, the critical crossover limit can be defined as R — >• oo, 
t, h — > 0, with t = t/G, h = t/Gh fixed. The variables that are kept fixed are the same, but a 
different mechanism is responsible for the change of the Ginzburg parameters: in short-range 
models we vary u keeping the range R fixed and finite, while here we keep the interaction 
strength u fixed and vary the range R. 

In this paper we will study a generalization of the model ( |1 . 1|) in the presence of medium- 
range interactions. We will show explicitly in perturbation theory the equivalence between the 
crossover functions computed starting from the continuum 4 model and the results obtained 
in the medium-range model. As a byproduct we will also compute the non-universal constants 
relating the two cases so that we can compare the field-theory predictions with the numerical 



results of Refs. |17 Hl9| , P^1 without any free parameters. The calculation will also give us analytic 
predictions for the large- R behaviour of the critical point. 

In numerical simulations (and also in experiments) the range R is always finite. It is therefore 
important to understand the behaviour of the corrections one should expect. We will show that 
for d > 2 the deviations from the universal behaviour are of order R~ d , provided one chooses 
the scale R as in Eq. ( |1.4| ). These corrections are non- universal and depend on all the details 
of the microscopic interaction: as a consequence they cannot be computed in the continuum 
field-theory framework. In two dimensions the behaviour of the corrections is not computable 
in the perturbative expansion around the mean-field solution. Indeed the perturbative limit we 
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consider — first we expand in 1/R at t fixed and positive and then we take the limit t — > 
- does not commute with the crossover limit. We conjecture that the corrections are of order 
log R 2 /R 2 as already indicated by the numerical work of Ref. |18|]. This behaviour has been 
explicitly checked in the large- N limit. 

The paper is organized as follows. In Sect. EJ we review the computation of the crossover 
functions in the field-theory framework, extending the calculations of Refs. In Sect. 

^] we introduce our model with medium-range interactions and in Sect. |] we show that the 
expansion around the mean-field solution is equivalent to the perturbative expansion of the 
4 field theory apart from non-universal computable renormalization constants. In Sec. [5] we 
discuss the corrections to the universal critical behaviour. In Sect. |6]we compare our theoretical 
results with the Monte Carlo data of Refs. 



17H19|, pEq] and we present a critical discussion of 



the crossover model of Refs. [T^,|23[. A detailed comparison with Monte Carlo results for the 
self-avoiding walk will appear in a separate paper [p4fl . App. [A] contains some details about the 
computation of integrals with medium-range propagators, while App. [B] discusses the medium- 
range model in the large- N limit verifying explicitly various results presented in the text. 
Preliminary results were presented in Ref. [EH . 



2 Critical crossover functions from field theory 



In this Section we report the computation of the various crossover functions in the continuum 
theory. As we described in the introduction, the idea is the following: consider the continuum 
4 theory 
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where <fi is an A^-dimensional vector, and consider the limit u 
t = t/G = tu~ 2 '^~^ fixed. In this limit we have 



0, t 
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(2.1) 
0, with 

(2.2) 
(2.3) 

(2.4) 
(2.5) 



are respectively the susceptibility and the (second- moment) correlation length. The functions 
F x (t) and F^(t) can be accurately computed by means of perturbative field-theory calcula- 
tions. There are essentially two methods: (a) the fixed-dimension expansion || [7|, ^(| , which is 
at present the most precise one since seven-loop series are available p7|P5[; (b) the so-called 
minimal renormalization without e-expansion [|ll|, [3(J which uses five-loop e-expansion re- 
sults |]3"T|, |3"2"f . In these two schemes the crossover functions are expressed in terms of various 
RG functions whose perturbative series can be resummed with high accuracy using standard 
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methods |33,33]. Here we will consider the first approach although essentially equivalent results 
can be obtained using the second method. Explicitly we have for F x (t) and F^(t): 



X exp 



— dx 



7(x) 



(D 2 exp 
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(2.6) 
(2.7) 



where t is related to the zero-momentum four-point renormalized coupling g by 



7(x) 



~ f ° I, dX u{x)W{x) 



exp 



dz 



1 



yo v{z)W{z) 
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j(x), v(x), and W(x) are the standard RG functions (see Refs. [P?|,PS[| for the corresponding 
perturbative expressions), g* is the critical valuef] of g defined by W(g*) = 0, and x*, £*, an d 
Uq are normalization constants. 

The expressions ( |2.6j ), ( |2.7| ) and (|2.8| ) are valid for any dimension d < 4. The first two 
equations are always well defined, while Eq. ( |2.8j ) has been obtained with the additional 
hypothesis that the integral over x is convergent when the integration is extended up to g*. 
This hypothesis is verified when the system becomes critical at a finite value of j3 and shows a 
standard critical behaviour. Therefore Eq. ( |2.8| ) is always well defined for d > 2, and, in two 
dimensions, for N < 2. For N > 2, one can still define t by integrating up to an arbitrary point 
g$. For these values of N, t varies between — oo and +oo. 



We will normalize the coupling g as in Refs. [|27| [28|, p3| so that in the perturbative limit 
— ^ 0, t — > oo, we have 



y 2(47r) d / 2 3 V 2 / 



Xit 



(d-4)/2 



(2.9) 



This implies that for y — > we have t « (y / '\ d ) 2 /( d ~ 4 \ and (£*)% ~ X*h ~ 1- 

For future purposes we will be interested in computing the expansion of F x (t) for t — ► oo. 
In two dimensions we have 



N + 2 , f 24nt 

1 + ^g 



24?rt 



iV + 8 
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24vrf 



(2.10) 



where 

D 2 (N) 
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7(2) 



247T Jo x { v(x)W(x) 
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+ x + N + 8 



(2.11) 



1 A review of present estimates of g* can be found in Refs. 
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For N > 2, g* should be replaced in D2(N) with the arbitrary point go that has been used to 
define t. 

Analogously in three dimensions we have 



F x (t) 



1 + 



N + 2 N + 2, (A8nVt S 
log 



24 7 rt 1 /2 288vr 2 t 

2 
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27N 2 + 52N - 472 D 3 (N) „ r , /2 , ~, 
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where 
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(2.13) 



In the large- N limit we have simply 

D 3 (N) 



N 2 



(48vr) 



+ 0(N). 



(2.14) 



The expansions ( p.lOj ) and ( 2.12|) are nothing but the standard perturbative expansions. For 
generic values of d we have 

F x (t) = W a nt~ Ami , (2.15) 

^ n 

where A m f = (4 — d)/2. For d = 4 — 2/n additional logarithms appear in the expansion. The 



reason of this phenomenon is well known p8|-|4*0|1 : the critical crossover limit corresponds to the 
massless limit of the standard 4 model which is known to have logarithmic singularities for 
these values of the dimension. 

The functions F x (t) and F^(t) can be computed using the perturbative results of Refs. 
P7L^8], a Borel-Leroy transform that takes into account the large-order behaviour of the per- 
turbative series |jil"],|i2||, and a standard resummation technique |33|]. Explicit expressions can 
be found for iV = 1,2,3 and d — 3 in Refs. 0,[7|. Here we compute F x (t) and F^(t) for the 



Ising model in two dimensions using the four-loop results of Ref. ||27|| . In order to improve the 
precision of the results, by means of appropriate subtractions, we have forced the resummed 
expressions to have the correct asymptotic behaviour for t — > 0, i.e. F x (t) ~ t -7 , F^(t) ~ t~ 2u , 
with 7 = 7/4, v = 1. The resummed expressions are well fitted by 



F x it) 



2 log* 0.9705 0.3513 0.01712 0.001822 

1 + -T^F- + + —■ FT- + =5 + 



-i 3/16 



3vrf 

logt 0.7377 
1 + + — ~- + 
2%t t 



t 2 

0.1635 

t 2 



t 3 

0.00390 



t 



t 4 

0.000275 



1/4 



(2.16) 
(2.17) 
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The resummation errors on F x (t) (resp. F^(t)) are less than 0.1% for t <; 0.2 (resp. t <; 0.1), at 
most 3% (4%) for smaller values of t. The fitted expression introduces an additional uncertainty 
which is less than 0.3%. 

The constants D^N) are non-perturbative constants since they require the knowledge of 
the RG functions up to the critical point. By means of a Borel-Leroy transform, working as 
before, we obtain in three dimensions 

D 3 (0) = 0.002473(6), (2.18) 

D 3 (l) = 0.002391(6), (2.19) 

D 3 (2) = 0.002204(5), (2.20) 

D 3 (3) = 0.001920(3). (2.21) 

£> 2 (1) = -0.0524(2). (2.22) 

The results we have described above apply to the high-temperature phase of the model. For 
N = 1 the critical crossover can also be defined in the low-temperature phase || and the 
crossover functions can be computed in terms of (resummed) perturbative quantities. Using 
the perturbative results[]of Ref. 0, we can compute F x (t) in the low-temperature phase of the 
three-dimensional Ising model. For \t | > 10~ 3 the numerical results are well fitted by 

0.00019\ 02372 



In two dimensions we have 




x 



t 



! 0.0140674 0.0021871 0.000150048 4.82764 • lO" 6 5.78079 • 10" 8 
W 2 + t W 2 + T 2 W 2 



(2.23) 



The error due to the approximate form given above is at most 0.3%, while the resummation 
error is less than 0.1%. For |t| < 10 -3 we can use 

F x (t) = It]- 1 - 2372 (0.06125 + 1.3455t 1/2 + 1.51525?- 124.395? /2 ) , (2.24) 

with errors of order 0.2%. The resummation error varies approximately from 0.1% to 3%. 

In the low-temperature phase we can also study the magnetization. Using the results of 
Ref. IBI we have 



F M (t) = (a)u 



-1/2 



0.3241 0.02751 0.001247 0.0000128\ a0868 



For 1*1 > 10" 4 , the error on this function is at most of order 0.5%. 



2 Some perturbative series contained small errors, see footnote 27 of Ref. pq| . 
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We have chosen the expressions (|2.23|) and ( [2.25| ) so that they reproduce the exact large-t 
behaviour of the crossover functions: 



F M (t) 



1 

W\ 



1 



1 



1 1 



+ 0(r 1 \o g t)j , 

+ 0(t- 1 \ogt)). 



(2.26) 
(2.27) 



Notice that the leading correction to F x (t) is negative so that F x (t) is non-monotonic. Such 
a behaviour has also been observed numerically JT8[] and predicted analytically |20| in two 
dimensions. 



3 The models and the critical crossover limit 

In this Section we will study the critical crossover limit in the presence of medium-range inter- 
actions by performing a systematic expansion around the mean-field solution. In particular we 
will show how to compute the critical crossover functions as a perturbative expansion in powers 
of t( d-4 )/ 2 . A general discussion of the mean- field limit can be found e.g. in Refs. |34|,f43"fl. 

Generalizing the discussion of the introduction, we assume we have a family of couplings 
J p (x) that are defined on a (^-dimensional cubic lattice and that are parametrized by p. For 
each coupling J p (x), we define, in analogy with Eq. ( |1.4j) , the following quantities: 

V p = £j„(s), (3.1) 

X 

R 2 = o^E^W. ( 3 - 2 ) 
zavR x 

U p (q) = 1 - ij p (g), (3.3) 
v p 

where J p (q) is the Fourier transform of J p (x). In the following we assume that there is a one- 
to-one correspondence between R and p, so that we will interchangebly think of the various 
observables as functions either of p or of R. 

Notice that because of the definition fl3.2|) we have for q — > 



and 



U R (q) « R 2 q 2 . (3.4) 

We assume that: 

(i) Jr{x) is uniformly bounded in x and R, i.e. \Jr(x)\ < C, independently of x and R; 

(ii) Vr and R 2 are finite for finite values of p. For p — ► oo, R 2 — > +oo and Vr ~ R d ; 

(iii) the system is ferromagnetic, i.e. U.R(q) > for all q ^ in the first Brillouin zone; 
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(iv) Uji(q/R) has a finite limit for R — > oo at fixed q, i.e. Hji(q/R) — > 11(g). This assumption 
is equivalent to requiring that Jr(Rx) has a finite limit J^x) for _R — > oo at fixed x. 
Furthermore we assume that Joo{x) is infinitely different iable in x = so that the integral 
/ d d qq 2n (l — 11(g)) exists for any n. Notice the following asymptotic property of 11(g): 
because of Eq. (|3.4j) we have 11(g) « g 2 for g — ► 0. 



We will investigate the properties of a class of models that generalize the theory defined in Eq. 
fll.lf) . We will consider a Hamiltonian of the form 

H = -— E j r( x i ~ x a) <p xi ■ ¥x 2 + ^E h x ■ Vx, (3.5) 

X±,X2 x 

where (p x are iV-dimensional vectors, and a single-site measure d^{ip x ) which we will assume of 
the form 

dfi{<p x ) = e- y ^d N ip x , (3.6) 

where V((p x ) is an even function (often a polynomial) of <p x which is bounded from below and 
that satisfies V(x) ~ \x\ p , p > 2, for |x| — - > oo. The partition function is defined by 

Z[h\ = (n^)e^. (3.7) 

J X 

The Hamiltonian we defined in the introduction is a particular case of the general model we 
discuss here. Indeed consider a family of domains Dr and define 

t ( \ - / 1 if x e D R , , , 

Mx) = { otherwise, (3 ' 8) 

and 

K(^) = cpl + A - l) , (3.9) 

with A > 0. It is easy to see that Jr(x) satisfies all the assumptions, as long as Dr satisfies 
some simple requirements, see App. |A.1| . To derive the relation between the two models (for 
simplicity assume h x = 0) let us rewrite the partition function as 



z[h\ = /n^ex P (^E M% - y) fx ■ <Py - E [pI + Kvl - 1) 2 ] } • (3-10) 



Then we rescale the field 



((3NV R X 1/2 



and define 



4(1 -2A) 96A 

In terms of we obtain again Eq. ( |1 . 1| ) with the coupling defined by Eq. ( |1.3| ). As we shall 
see in the following, the large- R limit is well defined if f3 goes to zero as 1/Vr. Thus r and u 
remain finite as R — > oo so that the large- R limit of the Hamiltonian ( p..l|) and of the model 
defined in this Section are identical. 
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For A — > oo, the model with partition function (|3.10|) reduces to the so-called N- vector 
model which can be seen as a particular case in which 



dii(<p s ) = 5{ifi - l)cT yv (3-13) 

The coupling Jr(x) defined in Eq. (|3.5| ) is ambiguous for x = 0. Indeed one can define a 
new coupling and a new single-site measure 

J R (x) = Mx)-K6 X , Q , (3.14) 
dfi{<p x ) ee d^ x )e N ^'\ (3.15) 

without changing the results. This step, which at first sight may appear trivial, can be inter- 
preted as a mass renormalization needed, as we shall see, to have a scaling theory. It is also 
needed from a mathematical point of view: in order to make the following derivation math- 
ematically rigorous, we will require K to be such that the Fourier transform Jn{q) satisfies 
Jr(q) > for all values of q. 

Let us now discuss the critical limit of the model. If we consider the critical limit with R 



fixed, Eq. (3"7F) defines a generalized 0(iV)-symmetric model with short-range interactions. If 
d > 2, for each value of R there is a critical point| Pc,R\ f° r P ~^ Pc,r the susceptibility and the 
correlation length have the standard behaviour 

Xr(P) « A x (R)t^(l + B x (R)t A ), (3.16) 
e R W) « MR)r 2 »(l+Bt(R)t A ), (3.17) 

where t = (j3 c ,R — P)/Pc,r an d we have neglected additional subleading corrections. The expo- 
nents 7, v and A do not depend on R. On the other hand, the amplitudes are non-universalQ. 
For R — > oo, they behave as fi~7|,[P8l 

A X (R) w A™R 2d ^'^ A \ A^R) ~ AfR i{2 - dv)/{4 - d \ 

The critical point f3 c> R depends explicitly on R. For R — >• oo we have /3 C) _r ~ l/i? d with 
corrections that will be computed in the next section. 



Let us now define the critical crossover limit. In this case we consider the limit [17, 18 



R — > oo, t — > 0, with _R M /( 4 d h ee i fixed. We will show perturbatively in the next section that 

XR = R- 2d/(4 - d) XR(P) - fM (3.19) 
^ ee R-VM&tf) - / c (t), (3.20) 

where the functions and /^(t) are universal apart from an overall rescaling of t and a 

constant factor, in agreement with the argument presented in the introduction. 

3 In two dimensions a critical point exists only for N < 2. Theories with N > 3 are asymptotically free and 
become critical only in the limit (3 — > oo. 

4 However some ratios of correction-to-scaling amplitudes are universal. For instance B^(R) / B X (R) is uni- 
versal and therefore independent of R. 
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There exists an equivalent way to define the crossover limit which is due to Thouless [15]. 
Let /J^r be the expansion of /3 Cj r for R — > oo up to terms of order R~ 2d /^~ d ) /Vr, i.e. such 
that 

lim R^^p-k (P C , R - P^f) = b c , (3.21) 



R^oo 

with \b c \ < +00. Then introduce 

t-- 



(3.22) 



It is trivial to see that in the standard crossover limit t = t + b c . Therefore the crossover limit 
can be defined considering the limit R — > 00, (3 — > /3£g^ with t fixed. The crossover functions 
will be identical to the previous ones apart from a shift. Thouless' definition of critical crossover 
has an important advantage. It allows the definition of the critical crossover limit in models 
that do not have a critical point for finite values of R: indeed, even if /3 c ,r does not exist, one 
can define a quantity (3^^ and a variable t such that the limit R — ► 00 with t fixed exists^. 

Moreover, as we shall see, (3^^ is known analytically: therefore, in the analysis of Monte Carlo 
data, no computation of (3 c ,r is needed and a source of errors is eliminated. 



4 Mean-field perturbative expansion 
4.1 General framework 



The starting point of our expansion is the identity []44|, 45 1 — we use matrix notation and drop 
the subscript R from Jr(x) to simplify the notation — 



exp 



N(3 



ipjip 



(detNpJ 



-N/2 f d 



N, 



(2tt) n / 2 



exp < N 



1 



+ 



(4.1) 



where <fr is another iV-dimensional vector field. The second ingredient is the single-site integral 
that defines the function A(<p): 



(4.2) 



where z is a normalization factor ensuring A(0) = 0. If we choose the single-site measure given 
in Eq. ( |3.6| ), we obtain the explicit formula 

roo / 9 \ N / 2 - 1 

ze NA W = ^N/2 I dxxN - le - V{x)+ N^/2 Ijf^NxW). (4.3) 

Jo \Nx\cp\) 



Notice that Eq. 



has a regular expansion in powers of <j) 2 , giving finally 



(4.4) 



5 This is the case of two-dimensional models with N > 3, see the discussion in Sec. 4.2, and of one-dimensional 
models with N > 1. In the latter case we can take /3^ p) = l/(a 2 V R ) where d2 is defined in Sec. n.ll 
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The first few coefficients are given by 



a 2 
a i 

«6 



3N 



N + 2 

15iV 2 



Nf A -{N + 2)fl 



(N + 2)(N + A) 



where 



fk 



N 2 h ~ SN(N + A)f A f 2 + 2(N + 2)(N + A)f\ 



J °° dx x N+k-i e -v{x)+NpKx 2 /2 



(4.5) 
(4.6) 
(4.7) 

(4.8) 



The results for the iV-vector model are obtained setting = 1 in the previous formulae. 

The coefficients a 2k depend on the various parameters that appear in the single-site measure, 
and on (3K. In the following we will assume V(x) to be independent of R, although considering 
-R-dependent potentials does not introduce any significant change in the discussion as long as 
liniR^oo V(x) exists and is finite. Under this assumption a 2k depends on R only through the 
combination f3K. In the following we will always consider the limit R — > oo with f3K — > 0, and 
therefore we will introduce 

a 2k = \im Q a 2k . (4.9) 

We will discuss the generic case[] in which a 4 < 0. If one tunes the parameters appropriately 
one can obtain a 4 = and a different critical limit. We will not consider these cases here. For a 
discussion in the large- N limit, see App. [B. 1| . In our expansion around the mean-field solution 
we will need the expansion of a 2 in powers of (3K. Explicitly we have: 



('■2 



a 2 + f3K 
+(3 2 K 2 



N + 2_ _ 2 

(N + 2)(N + A) 



120N 2 



N + 2__ ^ 

a 6 H a A a 2 + ^2 



+ 0(f3 3 K 3 



(4.10) 



Using Eqs. (fP] ) and (|4~2]) we have 



Z[h] oc jY[d<j) x exp jiV 



1 



J2 A{(j) x + k 



(4.11) 



The correlation functions of the (^-fields are obtained by taking derivatives with respect to h x . 
Using the equations of motion, it is easy to relate them to correlations of the 0-fields. For 
instance, for h T = 0, we have 



(<P* ■ <Pv) 



1 



J 



l/Z 



(4.12) 



6 Formally the discussion we will present requires only 124 ^ 0. However for 0,4 > the expansion we obtain 
would correspond to a </> 4 model with negative coupling. We therefore expect that potentials V(x) such that 
S4 > correspond to non-critical models. This can be checked explicitly in the large-TV limit for a potential 
containing a <fi and a 6 coupling. In the large-iV limit one can check explicitly that S4 < is a necessary 
condition in order to obtain the critical crossover limit, see Sec. B.l. 
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where the expectation value (<p x ■ (p y ) is obtained using the Hamiltonian ( |3.5|) , while (4> w ■ <p z ) 
is computed in the model ( [4.11|) . 

Now let us consider the (formal) perturbative expansion of the theory ( |4.11| ) with h x = 0. 
It corresponds to a scalar model with propagator given by 



AO?) 



N 1 - a 2 [3J(q) 



(4.13) 



and vertices 4 , 6 , that can be read from the expansion of the function A((f>), see Eq. (|4 
Let us now define 

1 K 

t = 1 + — , 4.14) 

a 2 V R (3 Vr y ' 

and consider the limit (3 — > 0, R — > oo with K and t fixed. If t > this formal perturbative 
expansion defines an expansion in powers of R~ d . To prove this result let us first rewrite the 
propagator as 



m = w 



(3K 



+ 



1 - U(q) 



Ai + A 2 (g). 



1 + a 2 f3K a 2 l + a 2 f3K U(q) + t _ 
Now let us consider a generic /-loop graph. It has the generic form 

fn«*l ot / u0- d n a(%) m^rs (e^ I - 



(4.15) 



(4.16) 



i<j 



where Q{N) is an iV-dependent constant. Let us now expand A(q r ) in the graph using Eq. 
(|4.15| ). We obtain a sum of terms that can be represented as graphs in which each line is 
associated to a propagator A 2 (q); these graphs are obtained from the original one contracting 
the lines corresponding to Ai. A generic term has the form 



(Ai 



n 

Kj 



dq 



; (2vr)^ 



1 nA 2 (^)ii( 2 ^(E^ 

i<j i 



(4.17) 



corresponding to an m-loop subgraph. Now notice that, because of the assumptions we have 
made at the beginning, A 2 (q/R) converges, for R — > oo at fixed q, to a function A 2 (q) given by 



A 2 (ff) 



i i-n(g) 

a 2 N U(q)+t' 



(4.18) 



that is integrable for all positive t. Then change variables in Eq. ( |4.17| ), setting q^ = PijR, and 
then take the limit R —> oo in the integrand, keeping t fixed. We obtain 



(AiT 

R dm 



n 

Kj 



(4.19) 



where the integration is extended over H m . The integral is ^-independent, and, as long as 
t is positive, it is finite. Since Ai ~ R~ d , the leading contribution of this graph behaves as 
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j^-d(m+n) _ Obviously m + n > I [the equality corresponds to those cases in which all lines with 
Ai belong to one- loop tadpoles], thus proving that the leading contribution of an /-loop graph 
behaves as R~ dl . 

It is immediate to see that the perturbative expansion is not uniform as t — > 0. Indeed 
for d < 4 infrared divergences appear, so that the coefficients of the expansion diverge as 
t — > 0. If one expands these coefficients in the limit t —>■ and chooses K appropriately - 
this corresponds to a mass renormalization and fixes the expression of — one obtains a 

new series which is a perturbative expansion in powers of t^ -4 )/ 2 with additional logarithms if 
d — A — 2/n, nGE The resulting expression can then be interpreted as the expansion of the 
critical crossover functions for t — > oo, provided that the limit R — > oo at t fixed followed by 
t — > is identical to the crossover limit. It is important to stress that this commutativity is not 
an obvious fact and indeed it is not true at the level of the corrections to the universal behaviour. 
In the following we will also show that graphs containing 6-leg (or higher-order) vertices can be 
neglected in the critical crossover limit. In other words one can simply consider the 4 theory 
obtained from Eq. (4.11). This will explicitly give the relation between the crossover functions 



computed in the medium-range model and those obtained in the field-theory framework of Sect, 
i 

We will first discuss the two-dimensional case, in which all these features can be understood 
easily, then we will present the general case. 

4.2 Two-dimensional crossover limit 

We wish now to discuss the critical crossover limit in two dimensions using the perturbative 
expansion of the model Q4.ll ) . Let us consider the zero- momentum correlation function 



G (m,n) = £ ((^ k(0 2 k . . . . (4.20) 



X2,—,X m ,yi,—,yn 



which contains m insertions of 2 and n fields 0, and its one-particle irreducible counterpart 

p(m,n) 

Let us consider a generic /-loop graph contributing to r^ m ' n \ and, to begin with, let us 
suppose that it does not contain tadpoles. We will be interested in the crossover limit t — > 0, 
R — > oo with R 2 t fixed. After expanding A(q) = Ai + A^q), we will obtain, apart from 
numerical factors, an expression which is a sum of terms of the form ( 4.171) . We wish now to 



show that, if K increases with R at most logarithmically, we can disregard all terms containing 
Ai. In other words, we can simply substitute A(q) with A.2{q) in the original graph. To prove 
this result, we begin by taking the limit R — ► oo keeping t fixed and rewriting the integral of 
Eq. ( f4.17|) in the form of Eq. (|4.19 ). We should now consider the behaviour of this integral 



for t — > 0. Simple power counting indicates that the integral is infrared divergent in this limit. 
Let us suppose that the graph associated to Eq. ( |4.19| ) does not contain tadpoles. In this case, 
in order to compute the leading infrared contribution, we can substitute A 2 (g) with its small-g 
expansion l/(a 2 iV(g 2 +£)) and then rescale q 2 = tp 2 , obtaining 



(A i) n (R 2 t) m ~ Mint I dp^ \ 
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where m is the number of loops and Mj nt the number of internal lines. Since, by hypothesis, 
the graph associated to Eq. ( f4.21|) does not contain tadpoles, the integral is finite and thus the 
prefactor gives its behaviour in the crossover limit. If the integral in Eq. ( |4. 19| ) is associated 
to a graph with tadpoles, it can be written as the product of an integral associated to a graph 
without tadpoles, and therefore behaving as in Eq. (|4.21| ), and a power of the one- loop tadpole 
integral. Now, using Eq. ( A.56Q , we have 



N 



d 2 q 
(2nf 



R 2 



a 2 (l + a 2 (3K) 



hA~t) « -7^— logt+— . 

47ra2 a 2 



(4.22) 



Therefore neglecting logarithms all contributions behave as _R- 2 ( 2m -- M ™t)A™ in the crossover 
limit. Now, if K increases at most logarithmically, Ai behaves as 1/R 2 modulo logarithms. 
Therefore all terms of the form ( |4.19| ) behave as 



{log R 2 ) p 

ft2(2m-M int +n) ' 



(4.23) 



for some p. Now, if N int is the number of internal lines of the original /-loop graph, N int = 
Mi n t + n. Moreover if the original graph does not contain tadpoles m + n > I for n > 1. Thus, 
for n > 1, 

2m + n- M^t = 2(m + n) - N int > 21 - N int . (4.24) 

Therefore contributions with Ai decrease faster and can be neglected. Thus the original /-loop 
graph we started from can be computed replacing A(g) with A 2 (g) and scales as 



1 

B? 1 



t 



(R 2 t) 



l-Ni; 



R2(2l-N int ) 



(4.25) 



without any logarithm. 

We should now consider graphs with tadpoles. As we already discussed these contributions 
can be written as the product of an integral associated to a graph without tadpoles, and 



therefore behaving according to Eq. ( 4.25 ), and a power of the one-loop tadpole diagram. 
Now, using Eq. ( f4.22p , we have 



N 



d 2 k 



A(k) 



(3K 



1 

1 + a 2 (3K + R 2 



1 i 1^° 2 
logt + — 

47ra 2 a 2 



+ o{R 



-2\ 



(4.26) 



The parameter K is free and we can take it i?-dependent at our will. We should also define 
PcR i i- e - the scaling behaviour of the temperature (3. To fix these two variables we require 



that the tadpole scales as R 2 for R 
achieved by taking 



oo without logarithms and that t ~ R 2 . This can be 



K 



log it 



c 



AttR 2 



V f 



R 2 



R 2 

N + 2a 4 1 
24ttN¥ 2 R 2 



R, 



\ogR 2 -a 2 (3V R 



(4.27) 
(4.28) 
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where Co is an arbitrary constant. In the derivation we have used the expansion of a 2 in powers 
of f3K, see Eq. (|4.10| ). With this choice we have 



— (t + c ) + o(R- 2 ), 



N 



d 2 k 



A(k) 



a 2 R 2 



-^log(t + c ) + C 2 -c 



+ o{R- 



where 



Co 



iV + 2a 4 
QN ¥ 2 



c - 



(4.29) 
(4.30) 

(4.31) 



Therefore, with this choice of K, in the limit R — > 00, (3 — > 0, with t fixed, the tadpole scales as 
R~ 2 , without logarithms of R 2 . It follows that all graphs scale as in Eq. ( |4.25| ) with additional 
logarithms of (tR 2 ). 

Now, a simple topological argument gives 



TX 1 

21 - N int = 2 - m - - + - - 4)\4, 

^ 2 k 

where Vk is the number of fc-leg vertices. Therefore, the graph scales as 

(R 2 t) l ~ N ™ 



(4.32) 



i? 4 



-2m-n 



(4.33) 



with additional powers of log(i? 2 t). In the limit R — > 00, t — > 0, with t fixed, the previous 
formula shows that graphs that have one or more vertices with six or more legs vanish faster 
for R —>■ 00 than graphs containing only 4-leg vertices. This means that in this limit we can 
simply ignore the </> 6 , 8 , . . ., terms in Eq. ( |4.11| ). In conclusion we obtain 



p(m,n) _ p(m,n)^>4— 2m— n (t -\- C ) 



2-m-n/2 



i 



(m,n) 



(4.34) 



where r| m ' ra ^ behaves as t~ l times logarithms of t. Therefore the loop expansion provides the 
expansion of the critical crossover functions in the limit t — > oo. For instance, if one considers 
the susceptibility, one obtains at two loops, 



f x (t) = lim X R 



a 2 



N + 2 a A 



R— >oo 



t + c 6A^ a 2 (t + c ) 



47T 



log(f + c ) + c - C 2 



(N + 2) 2 a 2 



36N 2 al(t + cof 



Air 



\og(t + c ) + c - C 2 



1 

Air V47T 



log(t + c ) +c -C 2 + 



2H 



N + 2 



O((t + c )- 4 ) 



(4.35) 



247r 2 ^'(3) 36- 



where H 

This result should not depend on Co- Expanding in t for t — > oo it is easy to check that no 
dependence remains and we can simply set Co = 0. 
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From the discussion we have presented, it is clear that the crossover functions can be ob- 
tained directly in the standard 4 theory with hamiltonian 



d d x 



2 + ^ 



(4.36) 



Indeed all graphs, except the tadpole, have been computed using the propagator of this theory. 
The tadpole has been dealt with differently as it should be expected: indeed this is the only 
ultraviolet divergent diagram. Therefore we have proved that we can rewrite 



(4.37) 



where F x is the crossover curve computed in the short-range theory and that was explicitly 
expressed in terms of RG functions in Sec. In field-theoretic terms fi x is the field renor- 
malization, s is related to the difference in the normalization of the fields and of the coupling 
constant and r is the additive shift due to the mass renormalization. Comparing the expansion 
(|2T0|) with Eq. J05|) we obtain 



T 



Nal 

a,i 
Na 2 2 

a<i 

n¥ 2 



N + 2 
24vr 



log 



2AnNa 2 2 



(JV + 8)|a 4 | 



D 2 {N) + 



N + 8 N + 2 



24vr 



6 



Co 



(4.38) 
(4.39) 
(4.40) 



where D 2 (N) is defined in Eq. ( j2.11| ). 
For the A^-vector model we obtain 



6 



N 



— log 

4?r & 



'4^(^ + 2)" 



+ C 2 



6D 2 (N) 1 N + 8 
N + 2 4^ N + 2' 



(4.41) 
(4.42) 
(4.43) 



As a final remark we wish to notice that we have followed here Thouless' approach to the critical 
crossover. For > 3 this is the only possibility since no critical point exists. For N < 2 the 
crossover functions defined in this way differ by a simple shift given by the constant r computed 
above. 

For N < 2 our results give also the large- R behaviour of the critical point. Since the critical 
theory corresponds to t — where t is the parameter appearing in the field-theory crossover 
functions, we have 



Pc,R 



1 



a 2 Vi 



R 



N + 2 a 4 1 , _ 2 r n _ 2 , 



(4.44) 



Notice that this result depends explicitly on the single-site measure through a 2 and a 4 , while it 
does not depend on the hopping coupling J{x). Indeed all the dependence on J(x) is encoded in 
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the expansion variable R 2 . For the N- vector model we obtain the simpler and iV-independent 
expression 

l 1 - 

(4.45) 



1 + ^w-^ R2 + w- +0 ^' 1] 



The presence of a logarithmic factor in (3 Cj r was already predicted by Thouless in a modified 
Ising model. 

4.3 Three-dimensional crossover limit 

The ideas of the previous paragraphs can be generalized to any dimension d < 4. For generic 
values of d one works as follows: first one considers the graphs with I loops contributing to the 
one-particle irreducible two-point function, where I satisfies 

2d 

dl < — , (4.46) 

computing their contribution in the limit R — * oo, (3 — > with iR 2d /( 4_d ) fixed. Then one fixes 
the scaling behaviour of K in order to cancel all the terms that scale faster than R~ 2d /( 4: ~ d \ The 
expression of /S^r is obtained requiring tR 2d ^ 4 ~ d ^ to be constant for (3 — > and R — > oo. The 
perturbative expansion becomes an expansion in powers of ?! 4_d ) /2 with additional logarithms 
when 2/(4 — d) is an integer, i.e. for d — 4 — 2/n. The reason for the appearance of these 
singular terms is well known |38|-p0|: the critical crossover limit corresponds to the massless 



limit of the standard 4 theory which is known to have logarithmic singularities for these values 
of d. 

Let us now discuss in more detail the three-dimensional case. For d = 3 we should consider 
the one-loop and two- loop graphs in the expansion of the two-point function (<p a (j) b ). If r 1 - 0,2 ^ is 
the irreducible two-point function at zero external momentum, we have at two loops 

r (o,2) „ I N + 2 aT {N + 2? a 2 TT 

1 ~ 1 - K/V R ~ "6AT a4Tl " ^6A^ G4TlT2 



N + 2 2m (N + 2)(N + 4) m2 



<T 3 - V J a 6 T? } (4.47) 



where 



18N 2 4 o 120]V 2 



Tr = Nf^l), M8) 
T 2 = N 2 [-^-A(q) 2 , (4.49) 



(2tt) 3 

T 3 = j ^^M<l)Mk)M<l + k). (4.50) 

Now, using Eq. ( |A.46|) , for (3K ~ 0(R~ 3 ) — we will show in the following this is the correct 
asymptotic behaviour — we have, 

t = pK i 1 i m 

1 l + a 2 (3K^ a 2 (l + a 2 (3K) 1>RU 

= -(1 - a 2 (3K) (l 1>R - a 2 f3K) - -^—^ (tR 6 Y /2 + o(i?" 6 ). (4.51) 
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The estimate of T 2 is easy and we find, cf . Eq. ( |A.49| ) . 



To 



8Tra 



■1/2 



+ o(R°). 



(4.52) 



Finally, using Eqs. ( |A.46|) and ( |A.72| ) we obtain for T 3 — we will call the associated graph 
"two-loop watermelon" — 



a 3 9 R 6 



32tt 



logt + C 3 



o(R- 



(4.53) 



where C3 is a constant given in Eq. ( A.73Q . 

Now let us consider the terms that scale slower than R~ 6 for R — > 00 at tR 6 fixed. Using 
the previous results we have 



N + 2ai 
6N V 2 



a 2 pK)(I liR -a 2 (3K) 



N + 2 



192tt 2 A^ 2 akRS 



log/r 



(iV + 2) 2 2 



36N 2 



<TiT 2 . (4.54) 



At first we neglect the last term proportional to T\T 2 . If K has a finite limit for R — > 00 and 
i ~ R~ 6 in the same limit, using Eq. (|4.14|) , we obtain f3 = l/(a 2 Vfi;)(l + 0(R~ 3 )). We can 
then determine K by requiring the expression fl4.54|) to be of order R~ 6 . This gives 



K = V, 



R 



Jl ' K+ 32^iV^ l0gjR + & 



where Cq is arbitrary. The scaling behaviour of (3 is obtained requiring t ~ 0(R 
introduce a variable t related to (3 by 



P 



1 



a 2 V R 



1 



N + 2a 4 
6AT if 



R 



32n 2 Na 2 2 R fi 



log it 



i? 6 



(4.55) 
If we 

(4.56) 



and define the critical crossover limit as the limit (3 — > 0, R —>■ 00 with t fixed, we find that 
indeed t ~ i?~ 6 . More precisely, using Eqs. Q4.14|) and (|4.10|) , we have 



tR b 



t + c + ci + o(-R°) 



where 



Cl 



(N + 2)(N + 4)a 6 (N + 2) 2 a 2 N + 2a 4 



120N 2 



On 



18N 2 of 



6N a 2 . 



t 2 xuiv u, 2 

a is defined in Eq. ( |A.40| ), and c is related to c by Eq. (^4.31|) . Notice that for R 
converges to a constant, cf. Eq. (|A.39|) , and (3 ~ l/(a 2 VR) as it was claimed before. 
Using Eqs. ( |4.56| ) and ( |4.55[ ) we can rewrite T\ as 



a 4 1 



32ir 2 Nc$R? 



\ogR 2 + 0{R~ 



a>4 



3N 



T 3 + 0(R~ 



(4.57) 

(4.58) 
00, K 

(4.59) 



Let us now go back to Eq. ( |4.54|) . We should now deal with the term proportional to T{T 2 that 
was neglected in the previous treatment. We have 



r( o,2) = (N + 2f a\ 
9216tt 3 A^ 3 a 5 



tR G 
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~ 1,2 ^ogR 2 + 0{R% 



(4.60) 



Because of the presence of the logarithm, this term does not scale correctly. However it is 
of order 1/t 1 / 2 , and terms of this order appear also at three loops. We will now show that 
this contribution is canceled exactly by the contribution of the three-loop graph in which the 
tadpole has been replaced by the two-loop watermelon. This graph is associated to the integral 

T* = iV 5 / ^A(p) 2 A(g)A(r)A(p + q + r). (4.61) 

It can be rewritten as 

T 4 = T 3 T 2 + 

+iV5 / (t% £)3 A(p) 2 A(g)A(r)[A(p + g + r) - A(g + r)]. (4.62) 



(2tt) 3 (2tt) 3 (2tt) 3 

In the crossover limit the last term can be easily computed using the technique presented in 
App. |A.3| . It is easy to see that we can neglect the contributions due to Ai, and that we can 



replace 

Thus, neglecting terms o(R~ 6 ), we have 

(tR 6 y 1/2 f d 3 p d 3 q d 3 r 



T 4 = T 3 T 2 + 



a^R 6 J (2tt) 3 (2tt) 3 (2tt) 3 

1 



(p 2 + l) 2 (g 2 + l)(r 2 + l) 



(p + q + r) 2 + 1 (q + r) 2 + 1 



(4.64) 



Thus, if we keep only terms of order R 6 logi? 2 , we have simply T4 « T S T 2 . Then including the 
combinatorial and group factors we find the total contribution 

( Ar + 2 ) 2 - T 3 T ^ {N + 2f ^ (iV + 2) 2 

" 108iV 3 4 4 lOSiV 3 " 4 3 2 + ( } ~ 36iV 2 a ^ + °i R h ( 4 -65) 



where we have used Eq. (|4.59 ). Comparing with Eq. ([4.54 ) we see that the logarithms cancel 



as claimed at the beginning. 

This calculation illustrates the general mechanism in three dimensions. Consider a graph 
with only 4 vertices that does not have tadpoles or two-loop watermelons as subgraphs. In the 
crossover limit the contribution of this graph can be computed neglecting A x and substituting 
A 2 (g) with its small-g expression. Following the argument presented in two dimensions it is 
easy to see that an /-loop contribution scales exactly as t x ~ 1 ^ 2 . The same argument we have 
presented in two dimensions proves also that the contributions of graphs with 6 , s , . . . vertices 
are suppressed and can be neglected in the critical crossover limit. Finally consider a graph 
with only 4 vertices that has tadpoles or two-loop watermelons as subgraphs. This graph 
generates non-scaling terms involving powers of log R 2 . However the logarithmic contribution 
associated to each tadpole is exactly canceled, using Eq. ( |4.59| ), by the analogous term which is 
associated to the graph in which the tadpole is replaced by the two-loop watermelon. The sum 
of all contributions scales as R~ 6 without logarithms of R 2 . Our discussion proves therefore 
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that the perturbative expansion corresponds to an expansion in powers of 1/t 1 ^ 2 of the critical 
crossover functions with additional logarithms of t. Explicitly for the susceptibility we obtain 
at two loops 



{UT})- 1 = lim X ^R % = a 2 t + 



N + 2a 4Tl/2 



N + 2a 4 



R^oo 

N + 2 a 2 
' 18iV 2 if 



t l/2 + a 2 ci + ^r^— a 2 



32tt< 



logt + Q 



24tiN a 2 

(N + 2f a 2 
11527T 2 N 2 al 



+ 



6N a 2 

+ o(r 1/2 logt). 



(4.66) 



As expected, the contribution proportional to a§ at two- loops disappears in the crossover limit: 
again only the 4 - vertex is relevant. The result is also independent on cq. 

The discussion we have presented shows also that the critical crossover functions can be 
computed in the standard continuum 4 theory. Therefore, as we did in two dimensions, we 
can write 

f x (t) = fi x F x [s(t + r)} } (4.67) 

where F x is the crossover curve computed in the short-range theory and that was explicitly 
expressed in terms of RG functions in Sec. ||. Comparing the expansion Q2.12j ) with Eq. ( |4.66|) , 
we obtain 



N 2 a\ 
N 2 <4 

772 ' 



(4.68) 
(4.69) 



iV + 2 a 2 

=jlog 



288vr 2 ;V 2 a 4 



AMNa 2 
(N + 8)\a 4 \ 



a 2 D 3 (N) a 2 9N 2 - 20N - 544 



+ 



TT'l 



N 2 



+ 



20736tt 2 A^ 2 



N + 2a A „ 
+ ci + — — — -a 2 



N + 2a\ 



C, 



18iV 2 a 4 3 ' 



6N a 2 ion u, 2 

where D 3 (N) is defined in Eq. Q2.13| ). For the N- vector model these equations become 
(N + 2) 2 



Vx 



36 
{N + 2) 2 



36 



1 



ln 2 (N + 2) 



log 



ln(N + 2Y 



36D 3 (N) 9N 2 - 20N - 544 



2C, 



(N + 8) J (N + 2) 2 576vr 2 (iV + 2) 2 N + 2' 



(4.70) 

(4.71) 
(4.72) 

(4.73) 



As remarked in the previous Section, Eq. (|4.56|) gives us the behaviour of the critical point c> r 
up to terms of order R~ 6 . We have 



a 2 V R 



N + 2a 4 
QN ¥ 2 



h,R + 



1 a 4 1 
32^ 2 n¥ 2 R® 



logR 2 



In the iV-vector case, this expression simplifies becoming 



Pc,R 



1 

V~R 



1 + h,R ~ 



16tt 2 N + 2R 6 



' log R 2 + 7 



R 6 



o{R- 



(4.74) 



(4.75) 
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The first correction to /3 Cj r was derived in Refs. while the presence of the logarithmic 

correction was predicted in a modified Ising model with long-range interactions by Thouless ]15] . 
Notice that for iV — > oo the logarithmic term disappears and that, using Eq. ( |2.14j ), r ~ 0(1/N) 
so that the resulting formula gives the exact large- N prediction for j3 c R . 



5 Corrections to the critical crossover functions 

In this Section we will study the corrections to the critical crossover functions. First we will 
present the results for the three-dimensional case, which is the easiest one, then we will dis- 
cuss the corrections in two dimensions. To derive the behaviour of these corrections we must 
introduce some additional hypothesis on the function Jr(x). We will thus assume: 

(v) E x (x 2 ) 2 J R (x) is finite. It follows 11(g) « q 2 + 0(g 4 ). 

(vi) In the limit R — > oo at fixed q, Un(q/R) has an expansion in powers of 1/R 2 . Explicitly 

OO 1 

(g/i?) = 11(g) + £ — II n (g). (5.1) 

n=l n 

Notice that, because of property (v), II n (g) ~ g 4 for q — > 0. 



5.1 Corrections in three dimensions 



We begin by discussing the three-dimensional case. We will show at two loops - - but we 
conjecture that this is true to all orders of perturbation theory — that the leading correction 
to the scaling behaviour is of order -R -3 , provided one appropriately defines t. In other words 
we will show that 

(5.2) 



X 



in the crossover limit for a suitable definition of t. This type of behaviour should be true for 
any dimension 2 < d < 4. It is obvious that an expansion of the form ( |5.2| ) cannot be valid 
generically. Indeed if (3^^ is such that Eq. ( |5.2j ) holds — t is defined in Eq. ( |3.22j ) — consider 
ttfi defined by 

(5.3) 



PNew,c,R — Pc,R \ L + /irL 



6— a 



with < a < 3. If tNew is the corresponding scaling variable we have iNew = t+AR~ a . Therefore 
the two definitions are identical in the critical crossover limit. However in the variable ^New the 
corrections are of order R~ a . 

Let us now go back to Eq. ( |4.47| ) , and again let us neglect at first the contribution propor- 
tional to T{r 2 . The expression for T\ appearing in Eq. (|4.51|) is valid up to terms of order R~ 9 
as it can be seen from the results of App. [A.2| . Using the expressions for K and j3, cf. Eqs. 



( |4.55| ) and ( f4.56| ), one finds that the leading correction in T\ is of order log R 2 /R 9 . Let us now 
consider T 3 . Using the results of App. |A.3| we find 



T 3 = -I 2 , R (t) + 0(R- 



■^logt + C 3 + § 



0(R- g ). 



(5.4) 
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Therefore with the definition of (3^^ appearing in Eq. ( [4.56 ) we would obtain corrections of 
order R~ 2 and R~ 3 logR 2 . We will now show that these corrections can be eliminated with a 
proper redefinition of the expressions of K and (3^^ . Considering for simplicity the N- vector 
case [we have a 2 = 1, and a 4 = —6N/ (N + 2)] we assume 
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(5.5) 
(5.6) 



where K±, K 2 , &i, and b 2 are constants to be determined. Then the corrections of order 1/R 8 
and log R 2 /R 9 in T\ and T 2 are the following: 
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where the dots indicate terms that scale as R~ 6 and i?~ 6 log-R 2 . To cancel the unwanted 
corrections we must require that the combination T\ — 2T 3 /(N + 2) is free of terms that scale 
as R~ 8 and _R _9 log-R 2 . In this way we determine the constants K\, K 2 , b\, and b 2 . Explicitly 



2F* 



N + 2' 



K 2 



2b 2 
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4vr 2 N + 2 ' 



(5.9) 



We must now consider the terms proportional to T\T 2 . As we already discussed before we must 
consider at the same time the diagram associated to T4. A simple analysis shows that Eq. 
( |4.64| ) has corrections of order 0(R~ 9 ) which are therefore negligible in the present discussion. 
Therefore, including the combinatorial and group factors we must show that 



T{T 2 - 



2T 2 T 3 
N + 2 



T 2 [T 1 - 



2T 3 
N + 2 



(5.10) 



is free of terms that scale as R~ 8 and R~~ 9 logR 2 . We have already shown that the term in 
parenthesis has this property. For T 2 , using Eq. ( |A.49| ), we can show that Eq. ( |4.52| ) is valid up 
to terms of order 0(R~ 3 ). Therefore the previous expression is free of the unwanted corrections. 
In conclusion we have proved at two loops Eq. ( |5.2j ). We conjecture this is true to all orders: 
graphs without tadpoles or insertions of the two-loop watermelon should have corrections of 
order R~ 9 , while terms of order R~ 8 and i? _9 logi? 2 that appear in graphs with tadpoles or 
two-loop watermelon insertions should cancel with the mechanism we presented above. At the 
order we are considering graphs containing vertices with more than four legs should still be 
negligible: at two loops the contribution proportional to a§ in Eq. flP7p scales as log 2 R 2 /R 12 . 
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Using the perturbative expansion we can compute the function g x (t) in the limit t — > oo. 
We have 

g x (t) = -f + 0(r^), (5.11) 



where E 3 is defined in Eq. ( |A.47| ) . Notice that this behaviour cannot be changed by modifying 



the definition of /3^ p ' ) . 

Let us now show that if one uses t defined using the exact /3 Cjj r one directly obtains the 
expansion ( |5.2|) . To prove this fact, assume the opposite and write 



X = f x (t) + ^h x (t) + o(R- a ), (5.12) 



with a < 3. For t — > and any value of R, x ~ ^ 7 - Therefore / x (i) ~ t~ 7 and h x {t) ~ t~ 7 in 
this limit. Now, it follows from our discussion that the term of order R~ a can be eliminated if 
we introduce a new variable t — t + AR~ a . Substituting in Eq. ( |5.12j ) and expanding in R~ a 
we have 

X = m + W J' X ® + ^h x (T) + o(R-). (5.13) 

Cancellation of the terms of order R~ a requires Af (t) + h x (t) = 0. However this relation 
cannot be true since f' x (t) ~ t _7_1 for t — > 0. Therefore h x (t) = 0. We have therefore showed 
that the variable t is a particularly good one, since it automatically eliminates a whole class of 
corrections to the leading behaviour. Another consequence of these results is that we can now 
estimate the order of the neglected terms in Eqs. ( |4.74| ), ( |4.75| ): the terms o(R~ 6 ) are of order 
R- s . 

In App. [B.l| we compute the function g x (t) for our general model in the large-iV limit. 
The graph of g x (t)/f x (t) for some particular cases is reported in Fig. [I]. We consider: (1) 
the A- vector model, (2) the potential V(<p) = N((p 4 — ip 2 ), and (3) the potential V((p) = 
N(ip G + ip 4 — ip 2 ). Although the function is not universal since it depends explicitly on various 
constants whose value is specific of the model one uses, the qualitative features are similar in 
all cases: g x (t)/f x (t) interpolates smoothly between the values for t = and t = oo. Notice 
however that the function is decreasing in the iV-vector model, while it is increasing in the 
other two cases. 

Finally let us notice that the result ( |5.2|) depends crucially on our use of R as scale and on 
the specific field normalizations used in the definition of our model. In general with an arbitrary 
scale p and arbitrarily-normalized fields we have 

X = A(p)f x (B(p)t) + y d g x (t) + . . . , (5.14) 
where A(oo) and B(oo) are non- vanishing constants. 
5.2 Corrections in two dimensions 

We wish now to discuss the corrections to the universal crossover curves in two dimensions. If we 
repeat the perturbative analysis we have performed in the previous Section we face immediately 
a difficulty. Working at one loop and using the results of App. |A.2| we find corrections to the 
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crossover functions of order log 2 R 2 /R 2 and log R 2 /R 2 . However, at variance with the three- 
dimensional case, only the former terms can be canceled with a redefinition of K and t: the 
terms proportional to log R 2 /R 2 cannot be canceled. Indeed let us suppose that 



K = ^ 



0= = TT 



Vr 

R 2 
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Vr 



±lo g R 2 + §lo g 2 R 2 + §lo g R 2 



(5.15) 
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At one loop we obtain 

p(o,2) = 



„ ,^log 2 i2 2 „ ,^\oz 2 R 



o(iT 4 ), 
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where the dots indicate terms that scale as R~ 2 . The coefficient T 2 (t) is given by 

1 



r 2 (t) = -6o 

This term can be canceled setting 

bo = 0, 

Let us now consider Ti(t). We have 
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(5.17) 



(5.18) 



(5.19) 



(5.20) 



where a\ and a 2 are determined by the low-momentum expansion of 11(g), cf. Eq. ( A. 53 ). This 
term does not depend on K\ or b\ and it is therefore impossible to eliminate it. Therefore at 
one loop we obtain correction terms of order log R 2 /R 4 . At two loops terms proportional to 
log R 2 /R 2 pop in and in general we have 



xf x {t)- 1 



R 2 



^og n R 2 g n {t). 



(5.21) 



n=0 



The presence of this infinite series of logarithms may indicate that perturbation theory does 
not provide us with the correct corrections and that a resummation of the perturbative series 
is needed. In other words the perturbative limit, R — > oo at t fixed followed by t — > may not 
commute with the crossover limit R — > oo, t —>■ at tR 2 fixed at the level of the corrections 
to the universal behaviour. This phenomenon is not new in two-dimensional models. Indeed a 



similar non-commutativity appears in the corrections to the finite-size scaling functions [^9jp0 



In the large- N limit, cf. App. p.2 , the corrections can be computed exactly and in this case 
one finds 



X 



/x(*) 



+ o(R- 



(5.22) 



However this simple behaviour may be due to the large- N limit. In general, as long as iV > 3, 
we do not expect a change in the exponent, but a more complicated behaviour of the logarithmic 
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corrections would not be surprising. By analogy with what has been found for the finite-size 
scaling corrections in Refs. [(49|,|5(|, we could have a behaviour of the form 



vf (fl-i^l | lOSR2 T 9 «® (5 23) 

Also for the Ising model it is unlikely that a new exponent appears. The numerical work of 
Refs. |jT7|,|I8| confirms this expectation: indeed they find that the corrections to scaling are well 
described in terms of a behaviour of the form ( |5.22|) . In these works A(t) and B(t) are assumed 
independent of t. Of course this is an approximation, but it is not surprising it works well, 
since these two functions should be slowly varying, as indicated by the large- N solution. 

One may wonder if the non-commutativity we have discussed above is peculiar of two- 
dimensional models. A simple analysis indicates that a similar problem should also appear in 
three dimensions if one considers the corrections of order R~ e . Indeed at this order T3 gives 
rise to terms log R 2 that cannot be eliminated by changing the scaling of K and (5. 



6 Discussion 

In this Section we wish to compare the analytic results obtained in the previous Sections with 
the numerical ones presented in Refs. |l7|- |T9i|22]| and discuss other approaches to the crossover 
problem. 



Let us first compare our results for {3 Cj r with the numerical determinations of Refs. [17,22 



These simulations are performed in the Ising model with coupling given in Eq. ( |3.8| ) and domain 
family 



D p = j:r:|> 2 <p 2 }. (6.1) 



Using the numerical results of App. [A] for the constants C2 and C3 and the numerical estimates 
of Sec. ^| for the non-perturbative constants D 2 (N) and D 3 (N), we obtain in two dimensions 
the asymptotic expression 

PcrVr « 1 + ^^gR 2 + °' 19 ^ 2 5(5) « 1 + -^(0.0796 log R 2 + 0.1975), (6.2) 

while in three dimensions 

Pc,rVr « 1 + 7,„ - — ^ log B? - 2^ZW , (6.3) 

where J^r is defined in Eq. (|A.37|) . Numerical estimates for selected values of p are reported 
in Table [l]. If one is interested in the expression of P c ,rVr up to terms of order o(R~ 3 ), one 
can replace Ii r with the asymptotic expression ( |A.39|) . In two dimensions Eq. ( |6.2| ) agrees 
approximately with the fit of the numerical data of Ref. [0. They quoteQ 

PcbVr w 1 + -^(0.076(3) log/? 2 + 0.172(7) ). (6.4) 



7 Notice the different normalization of R 2 : R 2 in Ref. Eal is four times our definition of R 2 . 



26 



To understand better the discrepancies we have considered 

Pc,R, approx 



A(R) = -l)R\ (6.5) 

\ Pc,R, exact / 



where /3 c ,i?,approx is the asymptotic form ( |6.2| ), while P c ,R,exact is the exact value determined 
in the Monte Carlo simulation. Asymptotically we should observe A(R) — > 0. However for 
the values of p used in the simulation A(R) shows a somewhat erratic behaviour. For p 2 = 
32,50,72,100,140, we have A(R) = -0.0812, 0.1649, 0.1660, -0.1943, -0.1124 with an error 
of approximately 5 • 10~ 4 due mainly to the uncertainty in Eq. ( |6.2|) (the error on /3 Ci # jexact 
is much smaller). Clearly these values of p are too small for the asymptotic expansion to be 
valid. Similar discrepancies are observed in three dimensions. The non-monotonic behaviour 
of the corrections appears to be a general phenomenon for the family of domains used in the 
simulations, and it is probably connected with the fact that the shape is not natural on a cubic 
lattice. Similar oscillation with R are observed in lattice integrals. For instance, from the 
results of Table [I] in App. |A.2| , one can see that the integral I^r does not have a monotonic 



behaviour even for p 2 « 10 3 . 

Let us now compare the results for the crossover curves. In Figs. 0, ^, f|, and [5|, we report 
the graph of the effective exponents 7 e fr, z/ e fj and /3 e ff defined by 

~ t df x (t) 

ues{t) = -2m-zr> (6 - 7) 

Affffl = (6-8) 
f M [t) dt 

for the Ising model in two and three dimensions, using the field-theory results presented in Sec. 
and the rescalings (|4.41 ) , ( 4.42j) , ( |4.71| ) , and ( |4.72j ) . In two dimensions we can compare our 



results for 7 e g(t) with the numerical ones of Ref. jTjJ. In Fig. || we report also the curve 

r\ - 3 1 + 0.339 1 1 / 2 ^ . 

7effW = 1 + 4 1- 0.115 tV2 + 4 .027r (6 - 9) 

whis is a rough interpolation of the numerical data. The agreement is very good, showing nicely 
the equivalence of medium-range and field-theory calculations. 

In Figs. | and || we report the results for 7 e g(?) and (3 e fi(t) in three dimensions. As already 
discussed in Ref. [|L9l] , in the high-temperature phase, 7 e fr(£) agrees nicely with the Monte Carlo 
data in the mean-field region while discrepancies appear in the neighbourhood of the Wilson- 
Fisher point. However, for t — > 0, only data with small values of p are present, so that the 
differences that are observed should be due to the corrections to the universal behaviour. The 
low-temperature phase shows a similar behaviour: good agreement in the mean-field region, and 
a difference near the Wilson-Fisher point where again only point with small p are available [5I| 



We can also compare the results for the magnetization. In Fig. ^| we report the combination 
2 — 7^ — 2/? e ff which should be compared with the analogous figure appearing in Ref. [fL9|| : the 
behaviour of the two curves is completely analogous. 
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We wish now to discuss a different approach to the crossover that has been developed in 
Refs. fl, 0, [13|, following the so-called RG matching ]53|, [54]] and that has been applied 
successfully to many different experimental situations [p|-[5|,p^,[52],[55|. These papers consider 
phenomenological parametrizations which are able to describe the crossover even outside the 
universal critical regime. Let us now introduce this model in the formulation of Ref . [J23] , which 
is intended to apply directly to our class of Hamiltonians. If t is the reduced temperature, one 
introduces two functions /t(t) and Y{t) defined by the set of equations 



(l-u)Y(t) 



c t tY(tf^ A , 



it 



At), 



1/2 



Y(t) 



u/A 



(6.10) 
(6.11) 



Notice that, although three non-universal constants q, A and u appear in these equations, Y(t) 
and K 2 (t)/ct depend only on u and on the combination Jctj A. The susceptibility is given by 



X 



clct 



t + 1 



Y{t) 



(7-l)/A, 



1 + 2/), 



(6.12) 



where 
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(6.13) 



where u* is a numerical constant, u* = 0.472, and c p is another normalization non-universal 
parameter. Notice that x _1 /( c p c t) depends only on u and ^/q/A, so that A or c t could be 
fixed to any value without loss of generality. In order to interpret the Monte Carlo results of 
Ref. Ref. |23j further assumes that c t and u scale as 



_ cto 



u 



R 4 ' 



(6.14) 



tR fixed, we obtain 
1 - 



u v 



Yo(t) 



2 1 + (2A - l)Yo(t) 



(6.15) 



In order to have the correct scaling of x, one should also set c p = c p qR. Then in the critical 
crossover limit t — > 0, R — > oo, with t 

r 1 = c 2 p0 c t0 tY (tf-^ A 

where Y Q (t) satisfies the equation 

1 - Y (t) = 

with 



Yo(t) 



1/2A 



(6.16) 



UqA 



(6.17) 



Eq. (|6.15|) defines the universal crossover function in this approach, the model-dependence 
being included in the constants a and ci c t o. In order to understand the accuracy of this 
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approach we can compare x obtained from Eq. Q6.15 ) with the very precise results of Bagnuls 
and Bervillier 0,0]. First of all let us compare the asymptotic behaviour for t — > and t — > oo. 
In the mean-field limit t — > oo we have 



X~ 



c p0 c t0 1 



92 



+ o(r 1 ) 



where #2 ~ 0.311. Using the results^ of the fit of Ref. 
we have 



u 



X 



C p0 CtO t 



(6.18) 

1.22/36, Cm = 1.72/6, A = tt, 

(6.19) 



with a ~ 0.062, to be compared with the exact result, cf. Eq. ( |4.66[ ), a = 1/(Att) « 0.0796. 
Notice that if we wish to reproduce the correct behaviour for t — > 00, we should also require 
c po c io — 1- Analogously for t — > we have 



X 



c P o c *o 1 + 



U V 



\ ofli-Vf (l - 9l a 2A t A 



0(t 



(6.20) 



where g± ~ 0.618. Using c 2 c t o = 1, we obtain numerically 



/T2A 



X" 1 = 2.49 1 7 (l - 3.42 t A + 0(t 

to be compared with 

XT 1 = (2.70 ± 0.04) f 1 (l - (4.0 ± 0.1) t A + 0(i 
obtained using the results of Ref. |7j and Eqs. ( |4.71| ), ( |4.72| ). Finally we report in Fig. |^ 

/x,phen(^) 



(6.21) 



(6.22) 



A*® 



/x,bb(*) 



(6.23) 



where / X)P hen(*) is given by Eq. ( p.l5| ) with the numerical values of Ref. El, and / Xi bb(?) is 
obtained using the expressions of Ref. and fixing the non-universal constants with the help 
of Eqs. ( J4.71[ ), ( |4.72| ) (therefore / x ,bb(?) does not have any free parameter). The agreement 
is overall good — the difference is less than 1.5% — except in a small neighbourhood of the 
Wilson-Fisher point where the difference increases to 8% as it can be seen comparing Eqs. 
( |6.21| ) and ( |6.22j ). Notice however that the region where the discrepancies are large is outside 
the domain investigated in the Monte Carlo simulation of Ref. ||19|| . Similar discrepancies were 
already observed in Ref. [|T3]j . 

Let us now consider the corrections to the leading behaviour. If we use the expressions ( |6.14|) 
we find corrections of order R~ A , in contrast with the theoretical analysis we have presented. 
However there is a simple modification that gives the correct corrections and that does not 
change the leading behaviour we have discussed before. It is enough to assume that, for R — > 00, 



_ uo 
U ^Rf 



CpO- 



(6.24) 



^Notice that we use a different normalization for R 2 : our R 2 is 1/6 of R 2 used in Ref. |H 
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Notice that this scaling of u is more natural, since, as we discussed in the introduction and in 
Ref. [ 20|1 , any coupling constant u should scale as R~ 3 in the crossover limit. Observe also that if 



one wishes to keep the interpretation of n(t) as an inverse correlation length, then n(t) ~ 

i.e. is a correlation measured in units of the interaction range. Using these rescalings, 

we can write _ 

X = f x (t) + ^9x(t) + 0(R- 6 ), (6.25) 

where g x (t) depends only on a 2 t apart from a multiplicative constant. By means of an explicit 
computation we obtain 



9x(t) 2( 7 -l)y (t) 



f x (t) (2A - l)Y (t) + 1 

~ u*v ~ 

(2A-l)Y (t) + l + —(l-Y (t)) 



1 Au*vY (t)(l-Y (t)) 

[(2A - i)y (t) + 1] 2 • 



(6.26) 



For t — > oo we have 



while, for t —>■ 0, we have 



^1 - - -0.45, (6.27) 

- gi a 2 *t*. (6.28) 



^ , , 2 A, A 



fx® 

The behaviour g x (t) / f x (t) ~ i A for i — > is not what one should expect in general, see Fig. [I] 
for an example in the large- N limit, and it is related to our assumptions on c p and c t . If we 
include 1/R 3 corrections, i.e. assume 

c P = c p0 + ct = do + j^, (6.29) 

then g x (t)/f x (t) would tend to a non- vanishing constant for t — > 0. For t — > oo, we should 
compare Eq. ( |6.27|) with the exact result g x {t)/ f x {t) = —E^/u derived in Sec. |5~1 . 



A graph of g x {t)/f x (t) as a function of a 2 t is reported in Fig. §. It shows a behaviour 
analogous to that found in the large- N limit, and also the numerical size of the corrections is 
similar. 

It should be emphasized that the function g x (t) is non-universal and that it cannot be 
determined in continuum field theory. Therefore the expression (OS) cannot be justified and 



represents some natural — but nonetheless totally arbitrary — generalization of the field-theory 
results. For the model at hand it provides a reasonable qualitative approximation, but this is 
not true for any model one can consider. For instance, in the large- N limit, the ratio g x (t) / f x (t) 
can be either decreasing or increasing, see Fig. 0, depending on the Hamiltonian. On the other 
hand, the function in Eq. ( |6.26| ) is decreasing for any choice of the parameters. Therefore 
this phenomenological extension is not even guaranteed to be qualitatively correct. If one is 
interested in phenomenological interpolations that can describe the crossover even outside the 
universal regime, one could proceed in a more straightforward way, distinguishing clearly what 
can be predicted using the field-theory approach (the limiting universal curve) and what is 
introduced phenomenologically (the corrections to the universal behaviour). For instance one 
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could use the essentially exact f x (t) derived from perturbative field theory and any arbitrary 
reasonable definition for the corrections depending on some parameters that could be fitted to 
obtain the best agreement between data and model. In this way one could also obtain good 
phenomenological interpolations of the numerical (or experimental) data. 

Finally we wish to comment on the role played by the terms [1 + A 2 /n(t) 2 ] appearing in 
Eqs. Q6.ll ), (|6.13|) . In our large-i? expansion they can be simply replaced by A 2 /nit) 2 with 



corrections of order R 6 (of order R 8 with the original scalings). Therefore, these terms 



that were introduced in Refs. [|T^, |52j in order to improve the behaviour in the mean-field 
region, represents a way to introduce additional corrections of order R~ 6 . In the analysis of the 
numerical results of Ref. [|19j they play little role, since 



= Y2 tY ^ 1)/A w °- 17 ^W (2 ^ 1)/A (6-30) 

and t < 0.05, Y(t) < 1, 6R 2 <; 1. In practice the leading term and the first correction g x (t) 
already provide a good interpolation of the data of Ref. fll9 |. 
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A Integrals with medium-range propagators 

A.l Lattice propagators 

In this appendix we will compute the quantity 

E eik ' x i (A.l) 

for two choices of interaction domain D C Z d . For integer p we define 

D« = [x e Z d : < p for i : 1, . . . ,d] , (A.2) 

Df = {xEZ d : X><| <p|. (A.3) 



In order to compare with the numerical results of Refs. [|17H19|,|22| we will be also interested in 
the following family of domains 



(A.4) 



We will not be able to compute (|A.1| ) for this class of domains. However we will obtain some 



numerical results that will be used in the main text. 
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Let us compute ( |A.1| ) for . The computation is trivial and we obtain 
where L — p + 1/2. Correspondingly we find 

1 AT 2 — 1 

V P = (2p + l) d = (2L) d , R 2 = -p(p + 1) = (A.6) 

Let us now consider the second case. The computation is now much more involved. The 
result can be expressed in terms of the determinant of two <i-dimensional matrices. Define 

A _ f (cosfc,-)* -1 fori = l,...,d-l, j = 1, ...,d; , A ^ 

11 ~ \ fd(kj) ioTi = d,j = l,...,d\ 

Bij = (cosAy-)* -1 . (A.8) 

Then 

Q (2) _ p ik-x = 

detB 



«s - e = • ( a - 9 ) 



lent 2 ) 

The result depends on the function of a single variable fd(k) given by 

f d {k) = -2 cos | (sin £;) d ~ 2 cos ^JfeL + (A.10) 

where, as before, L — p+ 1/2. Explicitly in two and three dimensions we have 

k 

f 2 (k) = 2cos-cosA;L, (A. 11) 
k 

fs(k) = —2 cos — sin k sin kL. (A. 12) 
Expanding in powers of k it is possible to compute V p and R. In two dimensions we obtain 

V p = 2p 2 + 2p + l = ^(4L 2 + 1), (A.13) 

V P R 2 = I p (l + p )(l + p + p 2 ) = 1(4L 2 + 3)(4L 2 -1). (A.14) 
In three dimensions we have 

V p = I(2p+l)(2p 2 + 2p + 3) = ^(4L 2 + 5), (A.15) 

V P R 2 = lp( p+ l)(2p+l)(p 2 + p + 3) = ^L(4L 2 -1)(4L 2 + 11). (A.16) 
For large values of p one finds 

v > - R2 - w^nry (A - 17) 
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To prove Eq. ( |A.9|) , let us suppose that the result has the form 



E«M [e ,kiL + Pi,de 



(A.18) 



i=l 



where and Z?,^ depend on but not on L. This Ansatz is a natural generalization of the 
result that can be obtained in two and three dimensions by direct computation. Using the fact 
that 



nJ3(*i, ...,k d )= £ e^n^^h, . . . , fc^o 



ik d nr)(^) 



n=—p 



we obtain = (— l) a! and the following recursion relations: 

sin kj 



cos ki — cos kd 

d-l -y 

i=i 



cos fc — cos k n 



sm 



•l) d sin 



ki + fcfi) 



(A.19) 

(A.20) 
(A.21) 



where in Eq. (TAT20|) i = 1, . . . , d - 1. Using a Xil = (2 sin(A;/2))~ 1 , Eq. flOOD , and the obvious 
symmetry under permutation of the labels of the coordinates, we obtain, for 1 < z < d, 



(sin fc,- 



id-2 



A'; 



COS — . 

Ilj=i jyi(cos fcj - cos fcj) 2 



(A.22) 



We should now prove that this expression solves Eq. ( |A.21| ) which is the consistency condition 
of the Ansatz (|A.18 ). Assuming d even, we can rewrite Eq. ( |A.21|) as 



f sin h 



xd-2 



v 

i= i ]l?=ijyi(cosA;, - cosfej 



. 



(A.23) 



Let us now use the following result: given x\, . . . , x n , consider the n-dimensional matrix M y - = 
x\~ l . Then it is easy to see that (in the mathematical literature this determinant is known as 
Vandermonde determinant) 



detM 



(-1) 



n(n-l)/2 



n— 1 n 

n n 

2=1 j'=£-f-l 



Now define the matrix 



Cij = 



(cos kj) 1 ] 



for i = 1, . . . , d — 1; 
for i = d. 



(saikj) d 2 

Using Eq. ( A. 24 ), it is easy to convince oneself that Eq. dAT23|) can be written as 

detC 



detS 



0. 



(A.24) 



(A.25) 



(A.26) 



Since d is even, one can express (sinA^-) 2 as a sum of cos 21 kj, with < % < d— 1. Thus the 
last row of C is a linear combination of the previous rows and therefore det C — 0. When d 
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is odd the discussion is analogous. We have thus proved that the consistency condition (|A.21| ) 
is satisfied. Therefore the Ansatz (|A.18|) with given by Eq. ( |A.22j ) is the solution of the 
recurrence relation ( A.19| ) that uniquely defines Vt p 2 \{k). Using again Eq. ( A. 24 ) we obtain the 
result ( |OD . 

If J p (x) defined in Sec. ^| is given by Eq. (p.8|), then J p (q) = £l Pt d(q)- We wish now to 
prove that fl Pt d(q) satisfies the properties mentioned at the beginning of Sect. |3|. Property (i) 
is obvious, while properties (ii) and (iii) depend on D p : they are satisfied if V p ~ R d and if, for 
any x G Z d , there exist[] x±, . . . , Xd 6 D p such that x = J2i a i x i, a i £ Z. To check the fourth 
property, define v = kR and consider the limit of Vt p ^{y/R) at fixed v. An easy computation 
for the domains and gives 



v R 



n (v) (A.27) 



d sin m } - 



with 

Q (v) = n— . (A.28) 

r= i ^ 

n (v) = -d\Y: 1 f 2 _ ' - » ( A - 29 ) 

u = v lim — . (A. 30) 

R^oo R V ' 

Notice that in two dimensions there is a simple relation for £Lq{v) f° r the two domains and 
D^. Indeed 



where 



For the family of domains VLq{v) was computed in App. A of Ref. |T7| finding 



d _\ I \u\ 



-d/2 



«o(«) =T {2 + 1 )[y) Jd/2iH) ' (A ' 32) 
where 

Ui = sj2{d + 2) v { . (A.33) 

From these expressions it is easy to see that 11(g) = 1 — Qo(o) satisfies condition (iv). 

For and it is also easy to show explicitly properties (v) and (vi) of Sec. [5]. Indeed 
in the limit R — > oo at v fixed we have 

1 oo 1 

-n p4 (k) = £ (A.34) 

V P n=0 U 

without odd powers of 1/R. Indeed, from the explicit results we immediately see that Q p> d(k) is 
even under the transformations k — > —k and L — ► —L. Moreover L 2 is an analytic function of 

9 Notice that it is not sufficient that V p ~ R d to ensure property (iii). For instance consider in one dimension 
the set Dp — {x — 2n, n E Z, \n\ < p}. 
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R 2 . Therefore for v fixed, Q p ^{k) is even in R, proving Eq. ( A. 34 ). Since we used the explicit 



expressions we computed before this proof applies only to the two cases we have studied. 
However we conjecture this is a general property of every family of domains that is cubic 
invariant. 

A. 2 One-loop integrals 

Let us now consider, for d < 4, the following class of one-loop integrals: 

r f 2X _ f l-II fl (fc) 

Jl ' fi(m) = 7 (2^ (n fl (fc)+m2)2 » (A ' 36) 

where Iljj(g) is a function with the properties mentioned at the beginning of Sect. ^| and |5|. 
We will be interested in computing these integrals in the crossover limit that corresponds to 
R — > oo, m — * 0, with m 2 R 2d ^ 4 ~ d ^ fixed. The integral (|A.35|) exists for m = only for d > 2, 



while the second one is always infrared divergent. In order to analyze the asymptotic behaviour 
of these integrals we will distinguish three cases: (a) d > 2; (b) d < 2; (c) d = 2. 

A.2.1 Case (a): 2 < d < 4 

For d > 2 the integral ( |A.35| ) is well defined for m 2 — > and thus we begin by studying 



h, R = lim /i^K). (A.37) 

m 2 ^0 

We wish to compute its asymptotic behaviour for R —> oo. Defining p = /ci?, we rewrite 

1 /• rfV l-n fl (p/fl) 1 r d d P l-n(p) 

y (2vr) d ^(p/i?) + m 2 ^ R d J (27c) d IT(p) ' 1 J 

using property (iv). The last integral can be extended over all R d . Thus we obtain 

h,R « ^ (A.39) 

with 

7 (27r) d n(p) V ; 

If additionally we assume that n^(g) satisfies properties (v) and (vi) of Sec. [|, we can easily 
prove that I^r admits an expansion of the form 

-I oo 

7 «=* + (A.41) 
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p 












DW 


3 





042971778 





043960387 





041600702 


4 





043202728 





043921767 





041279504 


5 





043319601 





043713672 





041423800 


6 





043386811 





043664053 





041387988 


7 





043428975 





043574469 





041394901 


8 





043457153 





043547206 





041384933 


10 





043491295 





043486698 





041398502 


12 





043510406 





043451767 





041386965 


14 





043522170 





043429899 





041392394 


16 





043529921 





043415345 





041389669 


18 





043535297 





043405187 





041392740 


20 





043539178 





043397824 





041391612 



Table 1: Estimates of R 3 Ii, R for various values of p for the three domains introduced in the 
text. 



In three dimensions, for D^ 2 \ D^ 3 \ we have respectively^ 

0.0435562069 



h,R 
h,R 
h,R 



R 3 

0.04336529 



R 3 
0.04139 



+ 0(R~ 5 ) 
0(R- 5 ), 



R 3 



+ o{R~ 



(A.42) 



Estimates of Ii jR for various values of R are reported in Table [I]. 



Let us now go back to I\^ R {m ) and let us compute the leading correction depending on m 



in the crossover limit. We rewrite 



d d q 



1 - Il R (q) 



(27r) d U R (q)(U R (q)+m 2 Y 
Setting q = v/R and using property (iv) of n#(g) we find the leading term as R — ► oo 

d d v 1 - U(v) 



(A.43) 



h,R( m 



m 



(27r) d n(f)(n(w) +m 2 



(A.44) 



where the integration is extended over R d . If H R (q) satisfies properties (v) and (vi) the neglected 
terms are of order m 2 /R d+2 [in three dimensions they are of order R" 11 ]. Now, for d < 4, the 
last integral is infrared divergent for m —* 0. Thus the leading contribution in the limit m — ► 

10 Luijtcn |B2] has noticed that an approximate expression of a can be obtained using the numerical results 
of Ref. jg7)- He writes a w 4.461imH^oo R 3 /Vr. For the three domains one obtains 0.0379, 0.0374, 0.0337 
respectively; these estimates are not very far from the exact values. 
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is obtained replacing 11(f) with v 2 in the denominator and with 1 in the numerator. We have 
therefore 



m 2 f d d v 



1 



R d J (27r) d v 2 (v 2 + m 2 ) 
m 2 r d d v 1 - U(v) 



R d J (2n) 



)(Il(f ) + m 2 ) v 2 (v 2 + m 2 ) 



(A.45) 



The first integral can be computed exactly, while in the second one we can simply take the 
limit m — > 0. We have finally 



( i -0 r ^+ B 4 +o(B "*"* /<4 " J)) ' 



where 



E* = - 



d d % 



(2tt)< 



1 - U(v) 



.,2^2 



(A.46) 



(A.47) 



n(w) 2 

Numerically, in three dimensions, we have for the interaction (|3.8|) and the domains 
and 



E 3 


= 0.058391 


for L>W, 


E3 


= 0.058545 


for D^ 2 \ 


E 3 


= 0.0635 


for D^. 



The computation of J\ t R{m 2 ) is analogous. We obtain 



(A.48) 



(A.49) 



A.2.2 Case (b): d < 2 

Let us now consider Ii^{m 2 ) for d < 2. In this case the integral is infrared divergent. Following 
the previous steps, we have 



h,B.{rn 



Analogously 



1 r d d v l-n(f) 



R d J (27r) d U(v) + m 2 



0{R' d - 2 ) 



d d v 



1 



R d J (2ir) d v 2 + m 2 R d J (2tt 



d d v 



1 



1 - U(v) 
Il(f ) + m 2 v 2 + m 2 



(47r)- d/2 rfl 



d\ m d 2 1 
+ 



2 J R d R d J (2tt) 



d d v 



1 - n(«) 1 
n(t>) u 2 



+ o(ir 



J^m 2 )^ (4vr)-^rf2-^ md " 4 



2 



+ 0(m c! - 2 i?- d ). 



(A.50) 
(A.51) 
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A.2.3 Case (c): d = 2 



Let us now consider the case d = 2. Also in this case the integral is infrared divergent for 
m? — > 0. However we cannot proceed as in the case d < 2, since the subtracted integral in Eq. 
( |A.50|) is ultraviolet divergent. First we set q = v/R and expand Hr(v/R) obtaining 

d 2 v l-U(v) 1+m 2 f d 2 v YlAv) 



T l 2\ I [ d 2 v 1 



i? 4 



(2vr) 2 (U(v) +m 2 ) 2 



0(R~ 



(A.52) 



(2tt) 2 n(u) + m 2 

Since IIi(t>) ~ v A , we can simply take the limit m — > in the last term. Let us now deal with 
the first one. Because of the lattice symmetry and of property (v) we have, for v — > 0, 

U(v) = v 2 + a x {v 2 ) 2 + a 2 v 4 + o(v 4 ), 



(A.53) 



where v 4 



jk. Then we rewrite 



+ 



d 2 v 1 - U(v) 
(2tt) 2 U(v) + m 2 
d 2 v f e~ v2 



d 2 v 



U(v) 



(2ir) 2 1 ) + m 2 v 2 + m 2 



+ 



- w2 [ai(t; 2 ) 2 + a2^ 4 ] 
{v 2 + m 2 ) 2 



(27r) 2 1 v 2 + m 2 



2\2 



{v 2 + m 2 ) 



(A.54) 



The first integral can be expanded in powers of m 2 neglecting terms of order m 4 logm 2 , while 
the second one can be computed exactly. We obtain finally 

d 2 V 1 - Ufa) 1 o \ Tft 2 . 2 , . 

(0 s 2TT , 2 = -7 lo S m + ^e) - — logm 2 + lE - 1 

(2^)^ ll(t>) + m z 47r 47T 



m 



16tt 



-m 



(4ai + 3a 2 ) (2 log m 2 + 2 7b + 1) + 



(2vr) 2 



U(v) 



U(v) 



d 2 v 



U(v) 



n(v) 



(v 2 ) 2 



+ 



2e- v2 [a 1 (v 2 ) 2 + a 2 v i ] 



(v 2 ) 



2\3 



0(R~ 4 log R 2 



(A.55) 



We obtain finally, up to terms o(R 4 ), 

h >R (™ 2 ) ' 
where 



1 , 2 , ^2 
log m H 



m 



AnR 2 



R 2 8nR 2 



(4«i + 3a 2 + 2) logm 2 + 



m 2 G x 



% (A.56) 



c 2 = -2* + 

4tt 



d 2 ^ 

(2VT) 2 



i - n(v) 



(A.57) 



and Gi and G2 can be computed from Eqs. (|A.52 ) and (|A.55[ ). For the interaction defined in 
Eq. (|3.8|) and for the domains D^ 2 \ we have respectively 



C 2 = -0.04578786 for and . 

C 2 = -0.05045 for D®. (A.58) 

The equality of C 2 for the domains and follows from the identity ( |A.31| ). Analogously 

1 



J ^ = W7p + gSp<*» + *•» + 2)('°S '" 2 + 1) - § + o(R->) 



(A.59) 
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A. 3 The two-loop integral 

We wish now to discuss the two-loop integral 

d d q d d k (l~U R (q))(l-U R (k))(l-U R (q + k)) 



(A.60) 



(2vr) d {2Ti) d (U R {q) + m 2 ){U R {k) + m 2 ){U R {q + k) + m 2 ) ' 

in the crossover limit R — > oo, m — > with m 2 R 2d ^ 4: ~ d ^ fixed. 

Since we wish to compute the integral in the large- R limit, we can rescale the internal 
momenta and use properties (iv), (vi) in order to rewrite /2,i?(^ 2 ) i n the form 

1 r d d q d d k (1 - n(g))(l - n(A;))(l - IT(g + k)) 



R 2d J (2n) d (2vr) d (11(g) + m 2 )(U(k) + m 2 )(U(q + k) + m 2 ) 
d d q d d k n x (g)(l-n(A;))(l-n(g + A;)) 



(A.ei; 



R 2d + 2 J (27r) d (27r) d (n(g) + m 2 ) 2 (II(A;) + m 2 )(II(g + k) + m 2 ) ' 
where the integrals are extended over 1R 2d . The neglected terms are of order 0(R 2d+i ) and 



0{R 2d+l 



in 



2d-4^ 



■ l2.R{ m2 ) is infrared divergent for d < 3 and therefore we will distinguish three 



cases: (a) d > 3; (b) d < 3; (c) d = 3. 
A.3.1 Case (a): 3 < d < 4 



For d > 3 the integral is finite for m 2 — > 0. In analogy with the discussion of Sec. |A.2.1| we 
have 



+ 



d d g d d k 
R 2d J (2Tr) d (2ir) d 



,2d-l 



(g 2 + l)(fc 2 + l)((g + k) 2 + 1) q 2 k 2 (q + k) 2 



+ o(R 



-2d\ 



(A.62) 



where 



! -2,R 



d d q d d k (l-n fl (g))(l-n fl (AO)(l-n fl ( g + A;)) 



(2vr) d (27r) d n R (g)n H (fc)II fi (g + k) 

I /■ d d g d d A; (l-n(g))(l-n(fc))(l-n(g + fc)) 



+ 0(i?- 2d - 2 ). (A.63) 



i? 2d 7 (27r) d (27r) 
A.3.2 Case (b): d < 3 



n(g)n(/c)n(g + k) 



For d < 3 the integral is infrared divergent for m — ► 0. As we did in Sec. A. 2. 2 we have 



h,R{m 2 
1 



m 



2d— 6 



d d g 



1 



i? 2d 7 (27r) d (27r) d (g 2 + 1)(F + l)((g + A;) 2 + 1) 
d d q d d k [(l-n(g))(l-n(fc))(l-n(g + fc)) 1 



R 2d J (27r) d (27r) d 



n(g)n(/c)n(g + fc) 



q 2 k 2 (q + k) 2 



+ o(R 



-2d\ 



In two dimensions we have explicitly 



f ( m 2 ) = _JL^ ( _J_^(i/3) _ J_) + 0(R~ 4 ' 
2 ' RK ' m 2 R A V24tt 2 ^ v 1 ' 36/ 



(A.64) 
(A.65) 
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A.3.3 Case (c): d = 3 

For d = 3 the integral is infrared divergent for m 2 — » 0. However we cannot proceed as in 
the case d < 3 because the subtracted integral in Eq. ( |A.64j ) is ultraviolet divergent. We will 
obtain the asymptotic behaviour using the same method we used in two dimensions to deal 
with Ii fi(m 2 ). We write 



i 

R G ' 
d 3 q d 3 k 
R 6 J (2vr) 3 (2tt) 3 



R 8 



(l-n(g))(l-n(A;))(l-n(g + A;)) e 



-q 2 ~k 2 -(q+k) 2 



U{q)U{k)U{q + k) 



q 2 k 2 (q + k)' 



where 



/ 2 cont (^ 2 ) 
F 3 



d 3 q d 3 k e - q 2 -k 2 ~{ q +k) 2 
(2tt) 3 (2tt) 3 (q 2 + m 2 )(k 2 + m 2 )((q + k) 2 + m 2 ) 
d 3 q d 3 k ni(g)(l-n(ife))(l-n(g + Jfe)) 



(2tt) 3 (2tt) 3 



U(q) 2 U(k)U(q + k) 



(A.66) 

(A.67) 
(A.68) 



where all integrals are extended over 1R 6 and terms of order 0(R 10 ) have been neglected. In 
order to compute /2° nt ( m2 )> let us define 



P(x, m) 



d 3 p e 



-p 2 +ip-x 



(27r) 3 p 2 + m 2 

2 sinhmlxl + e~ m ' a: 'erf ( m — 



57T \x\ 



- e mN erf m + 



\x\ 

T 



where 



Then, for m — ► 0, 



erf (a;) 



7T JO 



r o Cont fm 2 



1 



327T 2 

Collecting all terms we have 



log(9m 2 ) + 2^ E + 4vr / dx 



dt e 



x 2 P s (x, 0) 



(A.69) 



(A.70) 



1 1 

647T 3 X + 1 



h,R( m2 ) 



+ 0(m). (A.71) 



(A.72) 



where Cq is defined as 



C, = 



16-7T 



+ 



d 3 g ci 3 /c 
(2tt) 3 (2tt) 3 



roo r 

(log 3 + 7^) + 4vr / cfe x 2 P 3 (x,0) 
Jo 

(l-n(g))(l-n(A;))(l-n( g + A;)) e 



1 1 



647T 3 X + 1 

-9 2 -fc 2 -(<?+fc) 2 



n(g)II(A;)n(g + fc) 



q 2 k 2 (q + k) 2 



(A.73) 
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We have computed the constant C 3 
We obtain 



for the three different domains introduced at the beginning. 



B Large- N limit 



c 3 = 


-0.0127 


for D®, 


c 3 = 


-0.0127 


for D^ 2 \ 


c 3 = 


-0.0129 


for D®. 



(A. 74) 



In this Appendix we will compute the crossover functions in the large- N limit. In App. [B.l| we 
will discuss our general model for d > 2, while in App. [B.2| we will consider the two dimensional 
case for the iV-vector model. 



B.l Crossover limit for 2 < d < 4 

In this Appendix we will study the model introduced in Sect. 3 in the large- N limit following 
the strategy of Refs. |59-£jS|. We write V((p) = NW(ip 2 ) and study the limit iV —> 00 with f3 
and W{x) fixed. The basic trick consists in rewriting 



e~ NW ^~ JdpdXexp 



~K<p 2 ~ P) ~ NW{p) 



(b.i; 



The saddle point is given by the equations 



W\p) = \[3V R (l + m 2 ), p(3V R = y^[1 + h, R (™% (B-2) 

where Ii,r(w 2 ) is defined in Eq. (|A.35 ), while the two-point function is given by 

(<Po<Px) = 4 / tt rVl 2 " ( B - 3 ) 
p J H. R {p) + m 2 

The critical point corresponds to m 2 = and therefore the critical values of p and (3 satisfy the 
equations 

W'(p c ) = ^(3 c V Rj p c (3 c V R = 1 + 7 1)fl . (B.4) 

The critical value p c is the solution of the equation 

l+7i,fl= 2p c W'(p c ). (B.5) 



We will not need to solve Eq. (p.5| ) explicitly, we will only assume W(x) to be such that a 
positive solution exists. Since I\ R ~ R~ d , for R — > 00 we can expand 



n=0 

where 



Pc = J2 PcnT 1)R , (B.6) 
n=0 

2p c0 W'( Pc0 ) = 1, (B.7) 
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and p c i, p C 2, ... can be computed iteratively in terms of the derivatives of W(p) computed at 
p = p C Q. Correspondingly we obtain the expansion of f3 c for R — > oo. We have 

^ VR -^ + TT^w 2 h ' R+ (i + 2AW h > R + ° {R ] (R8) 

where W n = W^ n \p c0 ). Let us now discuss the scaling behaviour. Using Eqs. ( |B.2| ) and ( |A.4(j| ), 
we have 

( P P - PcPc)V R = A d — -m 2 + (E d - a)^ + o(m 2 R- d ) (B.9) 
where A d = (4vr)- d / 2 r(l - d/2) and a is defined in Eq. (|A.40|) . Now let us introduce^ 

S = p-p c , t = ^f^ = p c0 V R ((3 c -(3). (B.10) 

Pc,MF 



Using the first gap equation and Eq. ( |B.8| ) we obtain 

1 



8 



2p c0 W"( Pc ) 
where 



9 cr m 2 



+ 0(t 2 ,m\tm 2 J tR- 2d ). (B.ir 



l + 2p 2 PU 2 
2p 2 W 2 : 



;b.i2) 



'2 

and we have assumed W2 7^ 0. Going back to Eq. ( |B.9| ) we obtain finally 

So ° + 1#) (m2 - t] - Ad ^ - {Ed ~ a) w = 0(tm2 ' tR ~ 2d > m2R ~ 2d) > (R13) 

where 

"~ = t 1 " • (R14) 

Notice that this equation is also valid in the iV-vector model with Sqo = 1, si >o0 = 0. Let us 
now consider the critical crossover limit. We introduce 

m 2 = m 2 R 2d '^- d \ t ee tR 2d ^ A - d \ (B.15) 



and expand 



™ 2 = + Jtd9m(t) + o(iT d ). (B.16) 



The universal crossover function f m (t) is given by 

Soo[f m (t) -t}~ A d f m {tf~ 2) ' 2 = 0. (B.17) 



11 The normalization of t is chosen so that the results can be directly compared with those of Sect. [| One 
could have defined t=(/3 c — P)//3 C , as in Sec. |[ This choice does not change the leading crossover curve, but 
changes the corrections. 
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This results shows clearly the universality of f m (t). Indeed all the model dependence is included 
in the constant Sqo that can be eliminated with a proper rescaling of f m (t) and t. 

It is important to notice that Eq. ( |B.17| ) has a solution for t — > only if Sqo > 0. Indeed 
we can rewrite the equation as 



(B.18) 



For t — > 0, f m (t) — > 0, and thus we can neglect the last term in the right-hand side. Since 
the left-hand side is positive — we are considering the high-temperature phase — we should 
have Aa/soo < 0. Since A d < 0, we obtain > 0. Notice that, requiring the stability of the 
free energy in the limit R — > oo one obtains the condition W2 > 0, which, using Eq. ( |B.12| ), 
gives again > 0. It is important to remark, as we shall show below, that this condition is 
equivalent to the requirement a 4 < that was introduced in Sec. pTT| 
In three dimensions Eq. ( P-lTp can be solved explicitly finding 



fm(t) 



1 + 64vr 2 s2 t - 1 



647r 2 s^ 

It is easy to compute the correction function g m (t), obtaining 



;b.i9) 



g m (t) = fm(t) \{E d -a)- Xf m (t) (d ~ iy2 



Soo ^ A d f m {t)^ ^ 



;b.2o) 



where A = s 1]00 A d / s^. 

Let us now compute the asymptotic behaviour of f m (t) and g m (t)- For t — > 00 we have 



fm{t) 

g m (t) 



1 + ^-4)/2 + 0( ^-4) ) 



Eh 



a 



i + o{# d - 



4)/2n 



;b.2i) 
;b.22) 



For t — > we have 



fm(t) 

g m (t) 



2fit l 



d-2 



+ 0(t 



T2Aa 



l+0(t l 



where 7 = 2/{d - 2), A = (4 - d)/(d - 2), fi = ( Soo /| 

We can also compute the crossover function for the susceptibility. Since 



we have 



fx® 
9x® 



X 



PcO 
fm(t) 
1 



PcO 



9m(t)f X (t) 2 ~ —/*(*)■ 



a 

$00 



(B.23) 
(B.24) 



(B.25) 

(B.26) 
(B.27) 
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While f x (t) is universal apart from normalization factors, the function g x (t) is model-dependent. 
For this class of models, we have a three-parameter family of functions g x {t) parametrized, for 
instance, by s liOQ , Soo, and E d — a. The fact that the g x (t) depends only on three parameters 
should be due to the large- iV limit. For general values of N we expect g x (t) to be a non-trivial 
functional of the Hamiltonian. In Fig. [I] we report g x (t) / f x (t) for three different models: the 
iV-vector model and the theories corresponding to the potentials 

WV) = <^ + (^ 2 -l) 2 , (B.28) 
W(ip 2 ) = ^ 2 + (^-l) 2 + (^ 2 ) 3 , (B.29) 

using in all cases the coupling ( |3.8| ) and the domain ( |A.2j) . Notice that although the curves are 
quantitatively different, the qualitative behaviour is similar. 

We wish now to compare the results presented above with the explicit calculations of Sec. 
P . The basic ingredient is the large- N expansion of the integral 

dxX N+k - 1 e~ NW ^ = - [°° d yy^-l e N(lo g y/2-W(y)) ( g 3Q) 

2 Jo 

The saddle point equation is 

Y y = w '(y)' (B.3i) 

which is exactly the equation defining p c0 , cf. Eq. (|B.7|) . The large- N expression of Eq. (|B.30|) 
is obtained defining y = p c0 + z/VN and expanding in powers of 1/yN. In this way we obtain 

a 2 = Pco, (B.32) 

a 4 = -^o, (B.33) 

a 6 = -^(l2 Pc0 Wl + lQpl W*-W 3 ). (B.34) 

S oo VV 2 

These expressions are obtained assuming W 2 > 0. As it can be seen from Eq. ( |B.33| ), this 
condition is equivalent to requiring a 4 < 0. Using the previous expressions, it is easy to verify 
all the formulae reported in Sect. |Oj . Notice that in the N- vector case = p c0 = 1. 

All the considerations we have presented above apply if W 2 > 0. If W 2 = 0, but W 3 ^ 0, 
the leading behaviour of S changes and, for R — >• oo, we have 



5 -> n 2 _ JXI (m 2 - t) = qR d {m 2 - t). (B.35) 



R 

2 P >iy 3 

Thus, for R — > oo, keeping only the leading terms, we have 



,d-2 



If we scale 



(q + p c0 ) (m 2 -t)- PcoA d ^- w 0. (B.36) 

i2 4d/(4-d) m 2 t = (B.37) 

we obtain again the same crossover scaling function. To interpret these results we should notice 
that W"(p c ) = corresponds to the tricritical point. Therefore, here we are considering theories 



m 
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that for any finite R have a standard critical point, while for R — > oo converge to the mean-field 
tricritical point which is also a Gaussian point. We find that the scaling crossover functions are 
unchanged, although the scaling variables are different. In the framework of the introduction, 
these theories correspond to models in which the bare coupling u scales also with R as 1/R , 
so that the Ginzburg parameter G becomes G = (_R~ 2d ) 2 /( d ~ 4 ) = ^-4d/(d-4)_ 

In order to have different crossover functions we must consider a family a potentials such 
that W"(p c ) = for all R, i.e. consider theories at the tricritical point for any value of R. With 
our field normalization it is impossible to realize such a case, unless W((p 2 ) depends explicitly 
on R, i.e. W((p 2 ) = W(ip 2 ,R). Assuming therefore d 2 W(p c , R)/dp 2 = 0, a simple computation 
gives 

5 2 *^—(m 2 -t), (B.38) 
PcoW 3 

where W n = d n W(p, oo)/dp n evaluated at p = p C Q. Then we have 

/ 2 jA 1/2 r]—1 

I m — t\ . m - 



PcoA d —— = 0. (B.39) 



RescalingQ 

~2 = R 2d/(3-d) m 2^ l= R 2d/(3-d)^ (R4Q) 

we obtain the equation for the crossover function 



(m 2 - t) - A d m d ~ 2 = 0. (B.4i; 



valid, of course, for d < 3. One can go further and define multicritical crossover functions. If 
W"(p e , R) = ... W( n \p c , R) = 0, then, by a rescaling 

m 2 = R 2nd/{2+2n - nd) m 2 , t = R 2nd/{2+2n-nd)^ (R42) 

we obtain for R — > oo 

s (n) (m 2 - t) 1/n - A d m d ~ 2 = 0, (B.43) 

for a suitable constant s^ n \ 

B.2 Crossover limit in two dimensions 

In this Appendix we will discuss the critical crossover limit in two dimensions but we will 
consider only the iV-vector model, since it already exhibits all the general features. 
The gap equation is given by 

PVr= WTVT 2 = T^i 1 + Mm 2 )). (B.44) 

J {2n) z IL R (q) + m 2 1 + m l 



12 These rescalings can be derived in the formalism presented in the introduction starting from a Hamiltonian 
with a tfi 6 coupling. 
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The asymptotic behaviour of /^(m 2 ) is reported in App. |A.2.3| . Now define t from 

1 , „o CbQ , „o t 



R 



1 + 



4nR 2 



l 0gR 2 + ^i 0gR 2 + ^L___ 



R A R 2 



(B.45) 



Here we have introduced two free parameters ao and a\ that represent the possible ambiguity 
in the definition of f^R- Then assume for R — > oo at t fixed that 



m 2 = -^f m (t) + ^g m (t,logR) + 0(R- 6 \ogR 2 ). 
Using the gap equation we obtain for the leading term 

f m (t) + — log f m (?) = t + C 2 . 

47T 



;b.46) 



(B.47) 



Eq. ( B.47| ) defines implicitly the crossover curve for the correlation length f^(t) = 1/ f m (t) and 
for the susceptibility f x (t) = It is easy to check that this equation has the correct limits. 

For t — > — oo we have the standard asymptotic-scaling behaviour 



-4n\t\+4TTC 2 



1 + 0{e'^ 



fm(t) 

while, for t — > +oo, we have 
This expansion agrees with the results presented in Sec. EO 



i--^io g f+S + o(r 2 

4nt t v 



log?) 



(B.48) 



^B.49) 



For the correction term we obtain 



g m (t, log R) ■■ 

47T/ m (t) 

l + 4vr/ m (t) 



47T/ m (t) 
1 + 47T/ m (t) 



— (4«i + 3a 2 )f m (t) - a 



57T 



logi? 2 



— {4a! + 3a 2 )f m (t) log f m {t) + ai -G 2 - (G?i - C 2 )/ m (t) 

07T 



(B.50) 



Notice that in general there is no choice of ao which allow to cancel the logarithmic term. On 
the other hand, if one chooses Hamiltonians such that a± = a 2 = 0, the logarithmic term 
cancels. This class of Hamiltonians are called Symanzik tree- level improved p9" |. 
For t — > — oo we have 



-4tt t 



g m (t, log R) « 4Tcf m (t) -a logR 2 + ai~G 2 + 0(e 
while for t — > +oo we have 

g m (t,logR) w ^(4 ai + 3a 2 )t log(i? 2 t) - (d - G 2 ) t + O(logt). 

87T 



:b.5i: 



(B.52) 



Notice that if one takes a = and a x = G 2 , the corrections are strongly reduced in the limit 
t — > — oo. On the other hand, in the mean-field limit, the corrections do not depend on the 
scaling of (3. 
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Finally we wish to compute g x (t, log R). Using Eq. ( |B.25|) we have 




(B.53) 



For t 



oo we have 



g x (t,\ogR) = ---(4a! + 3a 2 + 2) \ogR 2 + 1 - rT (4ai + 3a 2 ) log? + 



G\ — C*2 



(B.54) 



which agrees with the perturbative result (|5.20|) . 
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Figure 1: Ratio 9x{t) I fxtf) ^ n ^ ne large- -/V" limit for three different models: (1) iV-vector model; 
(2) potential W( V 2 ) = V 2 + (<p 2 - l) 2 ; (3) potential W(<p 2 ) = V 2 + (<^ 2 - f) 2 + (<p 2 ) 3 . In all 
cases J(x) is given in Eq. Q3.8| ) with domain ( |A.2| ). Soo is a constant defined in Eq. ( |B.12| ). 
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Figure 2: Effective susceptibility exponent as a function of t in the high-temperature phase of 
the two-dimensional Ising model. The dashed line represents our interpolation of the numerical 
results of Ref. fI5f. In the mean-field limit 7 e fr = 1, while for t — > 0, 7 e fr = 7/4. 
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Figure 3: Effective correlation-length exponent as a function of t for the high-temperature 
phase of the two-dimensional Ising model. In the mean-field limit u e g = 1/2, while for t — > 0, 

^cff = 1- 



52 




1Q~ 6 1Q~ 4 1Q~ 2 10° 10 2 10 



Figure 4: Effective susceptibility exponent as a function of t for the high- (7^) and low- (7 eff ) 
temperature phase of the three-dimensional Ising model. In the mean- field limit jeg = 1, while 
for \t\ -> 0, 7 cfr « 1.237. 
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Figure 5: Effective magnetization exponent as a function of t in the low-temperature phase 
of the three-dimensional Ising model. In the mean-field limit (3 e s = 1/2, while for \t\ — > 0, 
(3 cS w 0.327. 
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Figure 6: Combination 2 — 7^ — 2(3 cS as a function of t in the low-temperature phase of 
the three-dimensional Ising model. In the mean-field limit this combination vanishes, while for 
t — > it is equal to the specific- heat exponent a 0.109. 
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Figure 7: Ratio f x , p hen(t) / f x ,BB(t) as a function of t. 
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